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SECTION  I 


Aiin nt  1 vp  MijI  1 1 1 <>vnl  16  Khnr.  Speech  Coder 
1 • 6 Introduction 

Under  the  eight  month  Adaptive  Multilevel  16  Kbps  Speech 
Coder  contract,  (DCA100-76-C-0002 ) , GTE  Sylvania  developed  software 
for  a high  quality  16,000  bit/sec  speech  transmission  terminal 
having  an  audio  bandwidth  of  3000  Hz.  This  software  was  designed 
to  operate  in  a full  duplex  mode  on  two  GTE  Sylvania  Programmable 
Signal  Processors  (PSP)  already  owned  by  the  Defense  Communications 
Agency . 

As  a benchmark  to  evaluate  the  voice  quality  of  the  Adaptive 
Multilevel  Speech  Coder,  GTE  Sylvania  also  developed  full  duplex 
software  for  Continuous  Variable  Slope  Delta  Modulation  (CVSD). 

This  CVSD  software  enabled  DCA  to  modify  and  optimize  the  internal 
operating  parameters  of  CVSD  and  to  experiment  with  new  CVSD 
algorithms . 

The  development  of  these  16  Kbps  voice  coders  was  motivated 
by  the  expected  deployment  of  a 16  Kbps  secure  voice  network. 

A necessary  part  of  any  transmission  network  is  the  subscriber 
terminal  or  telephone.  For  digital  communication  network  the 
transmitter  terminal  must  convert  the  analog  speech  to  digital, 
encrypt  the  data,  receive  incoming  data,  decrypt  it  and  reproduce 
the  speech. 

Unfortunately,  16  Kbps  is  too  low  a data  rate  to  send  voice 
without  the  use  of  data  reduction  algorithms  and  these  algorithms 
generally  introduce  audible  distortion.  The  goal  of  these  studies 
was  to  improve  the  performance  existing  16  Kbps  CVSD  terminals 
and  investigate  the  potential  of  new  terminals  employing  the 
principles  of  Adaptive  Predictive  Coding  with  Adaptive  Quantization 
( APCQ ) . 


1 


2. 

The  significant  results  of  this  contract  are  threefold: 

1.  The  voice  quality  of  CVf>D  can  be  Improved  through 
relatively  minor  changes  in  the  internal  parameters  or  through 
the  use  of  new  CVSD  algorithms. 

2.  The  APCQ  coder  provides  voice  quality  that  is  noticeably 
superior  to  the  best  CVSD  algorithms,  and 

3.  The  APCQ  coder  can  perform  well  even  in  the  presence 

3 

of  channel  errors  exceeding  1 in  10  . 

These  findings  show  that  it  is  technically  feasible  to  improve 
voice  transmission  at  1 0 Kbps  over  that  obtained  by  using  the 
present  CVSD  terminals.  Using  the  first  result,  relatively  minor 
changes  in  the  CVSD  telephone  subscriber  unit  can  provide  voice 
transmission  with  less  distortion.  The  second  result,  indicates 
that  a new  APCQ  coder  in  the  subscriber  terminal  could  provide  voice 
transmission  with  almost  imperceptible  distortion.  Although  the 
circuitry  needed  to  implement  these  adapjfc^e  multilevel  coders  is 
considerably  more  complex  than  that  needed  for  CVSD,  studies  performed 
under  this  contract  indicate  that  largescale  integration  techniques 
can  minimize  this  difference  in  complexity. 

In  fact,  we  show  that  the  APCQ  coder  can  be  implemented  around 
an  LSI  microprocessor  using  approximately  $700  in  parts  and  about  4 2 
watts  of  . Moreover,  the  processing  speed  of  this  microprocessor 

is  sufficiently  fast  that  it  can  analyze  and  synthesize  25  msec  of 
speech  in  16  msec.  Thus  40$  of  the  time  the  LSI  microprocessor  can 
perform  other  functions  such  as  synchronization,  encryption  and  telephone 
line  control. 
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3. 


Of  course,  CVSD  will  always  be  less  costly  to  Implement  than 
our  design  Tor  the  APCQ  coder.  Nevertheless,  our  design  for 
Implementing  the  complete  subscriber  terminal  with  the  APCQ  coder- 
may  cost  no  more  than  the  present  subscriber  terminal  employing  CVSD. 

This  equality  In  cost  occurs  because  analog  to  digital  conversion  of 
voice  is  only  one  of  many  features  of  the  subscriber  terminal.  Telephone 
line  control,  encypherment , and  synchronization  can  significantly 
impact  the  cost  of  the  subscriber  terminal.  Our  design,  based  around 
a microprocessor , provides  a cost  effective  method  of  implementing 
these  features  in  software  by  sharing  microprocessor  hardware  with 
the  APCQ  coder.  The  present  CVSD  terminal  requires  separate  hardware 
for  each  feature,  thus  raising  its  total  complexity  and  cost. 

Because  of  the  improved  voice  quality  provided  by  the  Adaptive 
Multilevel  coders,  further  work  should  be  performed  to  simplify  the 
algorithms, to  improve  their  hardware  design, and  to  test  their 
performance  under  realistic  operating  environments  more  fully. 
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SECTION  2 


H. 

Software  Simulations 

2.0  Introduction 

This  study  consisted  of  two  parallel  efforts  for  developing 
full  duplex  voice  processing  systems.  The  simulation  studies 
investigated  new  techniques  for  coding  speech  and  reducing  the  data 
rate  and  the  real-time  software  development  implemented  the  most 
promising  techniques  and  several  CVSD  algorithms  on  the  GTE  Sylvania 
PSP  for  further  evaluation  and  demonstration.  This  real-time  software 
was  delivered  to  the  Defense  Communications  Agency  for  operation  on 
their  two  PSP's  built  under  a prior  contract  (DCAIOO-T^-C-OOSS ) with 
GTE  Sylvania. 

Specifically,  GTE  Sylvania 

. Developed  simulations  of  adaptive  predictive  coders 
with  multilevel  Quantizers ( APCQ)  at  16  Kbps  using  several 
forms  of  error  signal  quantizers. 

. Studied  the  performance  of  these  APCQ  systems  in  the 
presence  of  channel  errors 

. Wrote  full  duplex  APCQ  software  for  the  DCA  PSP's. 

and 

Wrote  full  duplex  software  for  three  forms  of  continuous 
variable  slope  delta  modulation  (CVSD) 

The  remainder  of  this  section  will  describe  the  operating 
principles  of  APCQ  and  the  error  signal  quantizers  and  will  provide 
the  results  of  the  FORTRAN  simulations.  The  following  section  will 
than  describe  the  operation  of  the  real-time  CVSD  coder  and  the 


real-time  APCQ  coder. 
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2 . ] Operating  Principles  of  the-  APCQ  Terminal 

Figure  2-1  illustrates  the  principles  of  a speech  transmission 

1,2,3 

based  on  adaptive  predictive  coding  with  adaptive  multilevel 

^,5 

quantization  (APCQ).  The  APCQ  technique  estimates  or  predicts  the 
present  input  samples  sn  from  past  history  of  the  waveform,  that  is 

s = P ( s ,,s 

n n-1*  n-2 

The  residual  signal  or1  error  s.  lgnal  e =s  - s>  , along  with  the 

r i r i n 

function  P provides  sufficient  information  for  the  receiver  to 
regenerate  the  input  precisely.  In  general,  however,  the  error  signal 
is  quantized  and  is  not  sent  exactly.  This  quantization  distorts  the 
output  speech  so  the  choice  of  an  appropriate  method  for  quantization 
is  important  to  good  quality  speech.  A properly  designed  predictor 
will  produce  an  error  signal  having  less  dynamic  range  and  smaller 
variance  than  the  input.  Thus  a quantizer  operating  on  the  error 
signal  needs  fewer  bits/sample  than  one  operating  on  the  input,  and 
this  significantly  lowers  the  data  rate. 

The  optimum  predictor  depends  upon  the  statistics  of  the  input 
and  thus  the  predictor  parameters  adapt  as  the  input  changes.  These 
parameters  are  chosen  to  minimize  the  mean-squared  error  en  between  the 
predicted  samples  sn  and  the  actual  samples  sn , over  an  analysis 
Interval  or  frame  of  length  T which  is  typically  between  20  and  25 
ms  long.  This  frame  length  is  often  called  the  band  length.  Once 
calculated,  the  predictor  parameters  remain  fixed  for  the  entire 
Interval  but  change  from  frame  to  frame  as  the  input  statistics 
change.  Thus,  during  each  analysis  fram^ , the  algorithms  determine 
the  predictor  parameters  from  the  Input  data  and  form  the  error  signal 
sequence  by  filtering  the  input  dir  tallv. 


I 
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For  speech,  the  predictor  function  P was  broken  Into  two 


separate  smaller  predictors  I'l  and  IV,  as  1 I I ust. rated  In  Figure 
The  form  of  PI  acknowledges  that  speech  Is  often  quas  1 -per lod  1 c with 

A 

period  M.  Consequently,  PI  estimates  speech  as  where  the 

pitch  gain  a indicates  that  there  is  either  a gain  change  from  period 
to  period  or  that  the  speech  is  not  perfectly  correlated  with  period  M. 

For  calculation  of  the  pitch  period  M we  chose  the  average 
magnitude  difference  function  (AMDF) (6) 

AM1)F(  j)  = Z I*.  - s.j  | 

».-l 

as  a computat ionally  efficient  method  for  its  estimation  pitch  M.  This 
function  has  no  multiplications  and  requires  little  numerical  scaling. 

It  also  has  a sharp  null  for  that  value  of  J which  usually  corresponds 
to  the  period  M. 

In  speech,  not  all  pitch  periods  are  possible  and  thus  the 
terminal  calculates  the  AMDF  for  those  values  of  J corresponding  to 
pitch  frequencies  between  70  and  3^0  Hz.  Accurate  pitch  extraction 
requires  additional  logic  to  test  for  false  or  multiple  nulls,  but  since 
later  processing  can  partially  correct  for  occasional  pitch  errors,  the 
period  M is  set  to  that  value  of  j corresponding  to  a minimum  value  of 
AMDF.  To  speed  up  the  processing  further,  the  AMDF  calculation  in  the 
APCQ  system  sums  over  every  third  value  of  the  1 ripu  t . This  causes 
spectral  aliasing  and  occasional  shifts  in  nulls  in  the  AMDF,  but 
listening  tests  indicate  no  impairment  in  the  quality  of  the  resulting 
speech.  Furthermore,  in  calculating  the  AMDF,  the  algorithms  form 
partial  summations  which  are  then  scaled  and  added  to  form  the  final 
summation.  Overflows  are  always  clamped  to  the  largest  number  of 
the  processor  but  this  does  not  affect  the  minimum  values  needed 
for  pitch  extraction. 
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To  minimize  the  total  squared  error  E = ^ - a3n_;vi)‘5i 

can  di  rferent,  1 ate  K with  respect  to  a and  sot  the  result  to  zero 


giving 


a = 52  *«*"  «/Z  *»*  '* 


where  T = the  frame  baud  length  In  samples. 

For  periodic  sounds,  such  as  vowels',  a Is  close  to  unity,  but  for 
noiselike  consonants  which  have  little  correlation  M samples  apart,  a 
is  near  zero. 


The  reduced  waveform. 


v.  = s„  — as.-.u, 


or  first  error  signal  still  contains  sufficient  redundancy  such  that 
a second  Nth-order  predictor  P2  can  further  reduce  the  output  signal 
energy,  especially  if  the  speech  is  not  periodic  or  if  the  period 
is  estimated  incorrectly.  This  second  predictor  uses  a weighted  sum  of 
N past  samples  of  the  speech  waveform  to  form  the  estimate 


i,  = £ a.fn-. 


where  the  a^’s  are  chosen  to  minimize  the  squared  error 

= Z (»,  - ».)> 


and  are  given  by  the  solution  to  the  matrix  equation 


<t>a  = C 


where 


and 


This  is  a symmetric  Toeplitz  matrix  because  the  elements  along 

any  diagonal  parallel  to  the  principal  diagonal  are  identical.  Efficient 

solutions  exist  that  supply  the  a.'s  and  the  mean-square  energy  U in 

7 1 

the  difference  signal.  In  addition,  because  the  elements  of  the  matrix 

arise  from  an  autocorrelation  function,  the  stationary  matrix  solution 

for  the  a1’s  will  yield  a stable  filter  during  synthesis,  with  the 

recursive  filters  shown  in  Figure  2-2.  Unfortunately,  the  a^'s  do  not 

make  good  transmission  parameters  because  quantization  or  errors  in 

transmission  can  cause  the  poles  of  the  receiver  filter  given  by 

1 / ( 1 — P 1)  (1-P2 ) to  move  outside  the  unit  circle  in  the  Z plane,  causing 

unstable  waveforms.  Consequently,  auxiliary  parameters 

called  partial  correlation  (PARCOR)  coefficients  K.  (which  are  the 

1 A 

negatives  of  reflection  coefficients  discussed  by  Aral)  are 
calculated  from  the  a^'s  and,  as  long  as  these  PARCOR  coefficients 

ZT  7 

have  magnitude  less  than  unity,  system  stability  is  assured.  ’ 
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rn  the  actual  terminal.  Instead  of  quantizing  the  error  signal 
as  Figure  2-2  depicts,  the  analyzer  is  reformulated  and  the  quantizer 
placed  inside  the  filtering  loop  as  shown  in  Figure  2-3.  Without  this 
quantizer,  this  configuration  has  the  transfer  function  (1-P1)(1-P2) 
as  discussed  earlier.  With  this  new  formulation,  however,  we  eliminate 
the  effects  of  small  errors  in  accumulation  caused  by  the  error 
quantizer.  This  occurs  because  the  only  input  to  the  predictor  is  the 
quantized  error  signal  and,  consequently,  the  predictor  at  the  analyzer 
generates  the  same  predicted  sequence  as  the  receiver  does  in  the 
absence  of  transmission  errors.  Since  it  compares  this  predicted 
signal  with  the  original  input,  the  analyzer  can  correct  for  distortions 
caused  by  the  error  signal  quantizer,  altering  the  error  signal 
sequence.  This  would  not  be  possible  without  the  feedback  arrangement. 

There  are  numerous  methods  for  quantizing  the  error  signal, 
and  Jayant  provides  a good  discussion  of  adaptive  techniques  which 
allow  at  least  four  levels  or  two  bits  per  error  sample  (12).  A major 
portion  of  this  study  effort  was  devoted  to  investigating  various 
quantizers  and  these  will  be  discussed  in  Section  2.2  . 

After  computing  the  pitch  gain  a,  the  period  M,  and  the  PARCOR 
coefficient  , the  analyzer  digitally  filters  the  speech  and  quantizes 
the  error  signal.  Then  the  K^'s,  a,  M,  the  quantized  error  sequence 
en,  and  a normalization  factor  NORM,  needed  to  perform  the  operations 
in  fixed-point  arithmetic,  are  quantized  and  sent  to  the  receiver  where 
the  predictor  coefficients  are  regenerated  in  an  iterative  fashion 
by  computing 

<”  = - a, 

« K, 

where  represents  a,  on  the  ith  iteration. 


./'  = i. - |, 
i = i. 
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The  synthesizer  then  creates  an  output  time  waveform 
that  Is  both  Intelligible  and  pleasing  to  listen  to.  Moreover,  this 
output  ! ::  reasonably  I nsenr.  1 1 1 ve  to  errors  1 r > pitch  extraction  because 
the  1’?  predictor-  on  the  reduced  waveform  arid  the  quantization  of  the 
error  signal  can  partially  compensate  for  wrong  pitch  values  used  in 
the  first  predictor  PI.  In  fact,  if  the  pitch  period  M incorrectly 
doubles,  as  often  happens  in  practice,  then  predictions  made  by  the 
first  predictor  PI  are  made  from  two  periods  before  instead  of  the 
previous  period.  This,  is  riot  a serious  error.  If'  other  incorrect 
values  for  M are  chosen  then  different  values  of  u arid  PARCOR 
coefficients  are  also  computed  to  compensate,  in  part,  for  this 
error.  Finally,  the  error  signal,  even  though  coarsely  quantized, 
carries  considerable  information  about  the  true  pitch,  should  this 
pitch  be  incorrectly  measured.  Thus  the  APCQ  speech  process  can 
function  without  severe  degradation  in  noisy  acoustic  environments  and 
on  many  different  speakers  where  accurate  pitch  extraction  is  difficult. 

The  APCQ  process  can  be  modified  to  obtain  other  forms  of 
speech  coders  currently  being  studied  and  discussed  in  the  literature. 
These  modification  can  yield  Adaptive  Pulse  Code  Modulation  Coders 
(APCM),  Adaptive  Differential  PCM  Coders  (ADPCM),  and  Adaptive 
Predictive  Coders  (APC).  These  configurations  for  a 16  Kb/s  voice 
coder  are  shown  in  Figure  2-^1.  The  major  differences  between  the 
different  forms  of  coders  are  in  the  complexity  of  the  predictor  and 
in  the  form  of  the  error  signal  quantizer.  For  example,  in  Adaptive 
PCM  (APCM),  there  is  no  prediction,  only  the  error  signal  quantizer 
adapts  to  the  input  speech.  In  Adaptive  Differential  PCM  (ADPCM) 
the  predictor  uses  a weighted  value  of  the  previous  sample  to  form 
its  estimate.  The  weight  never  changes  and  its  value  is  typically 


near  0.9. 
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There  In  no  pitch  prediction.  Again,  an  adaptive  quantizer  is 
incorporated  which  adapts  to  the  input.  As  already  discussed,  the 
adaptive  predictive  coder  with  adaptive  quantization  (APCQ)  has  two 
fairly  sophisticated  predictions:  the  coefficient  predictor  with 

typically  A weights  and  a pitch  predictor  with  pilch  gain  u.  The 
pitch  loop  can  be  easily  eliminated  without  reconfiguring  the  coder 
by  setting  the  pitch  gain  a to  zero.  The  Adaptive  Predictive  Coder 
(APC)  is  identical  to  the  APCQ  coder  except  for  the  form  of  error 
signal  quantizer  employed.  Its  quantizer  has  only  two  levels  whose 
amplitude  changes  once  per  frame  interval.  Its  data  rate  is  typically 
6400  bps.  The  APCQ  coder  employs  an  error  signal  quantizer  having 
3 or  more  levels  whose  amplitudes  change  either  once/frame  of  after 
each  sample.  The  data  rate  is  correspondingly  higher,  typically 
9600  to  16,000  bps. 


2.2 
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Error  Signal  Quantizers 

The  design  of  the  error  signal  quantizer  significantly  affects 


the  voice  quality  of  the  synthesized  speech.  Under  this  contract 
we  have  investigated  several  types  of  quantizers  but  considerably 

1 

more  work  needs  to  be  performed  in  this  area  to  find  a simple 
quantizer  which  performs  well  with  channel  errors  and  which  will  permit 
us  to  reduce  the  complexity  of  the  adaptive  predictor. 

The  significant  feature  of  all  the  quantizers  investigated  is 
the  adaptively  changing  size,  Q.  For  example.  Figure  2-5 
illustrates  an  8 level  quantizer.  This  quantizer  takes  an  input  signal 
and  converts  it  to  one  of  eight  discrete  levels.  The  number  of  this 
level  (000  to  111  in  binary  notation)  is  sent  to  the  receiver  which 
converts  it  back  to  the  desired  amplitude  value.  Since  the  curve 
describing  the  input/output  characteristics  of  the  quantizer  looks 
like  a staircase  of  step  high  Q (for  uniform  quantization),  Q is 
commonly  referred  to  as  the  step  size.  Until  recently,  these 
quantizers  were  memoryless,  and  once  calculated,  Q never  changed  even 
if  the  signal  statistics  changed.  Unfortunately,  if  the  input  to  the 
quantizer  becomes  large,  the  quantizer  tends  to  saturate  and  not 
describe  the  input  well.  In  a similar  fashion,  a very  low  level  input 
may  force  the  quantizer  to  always  output  its  minimum  step  size  and  not 
track  the  small  variations  in  the  input.  What  is  needed  is  a quantizer 
having  memory  which  will  increase  the  step  size  for  large  signals  or 
signals  with  large  variance  and  decrease  it  when  the  input  signals 
become  small  or  the  variance  becomes  small. 

All  the  error  signal  quantizers  studied  in  this  contract 
changed  their  step  sizes  at  regular  intervals.  By  changing  the  step 


. 
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Figure  2-5.  8-Level  Adaptive  Quantizer  (Jayant  Quantizer) 
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sizes  in  this  manner,  the  quantizer  adapts  to  the  changing  variance 
of  its  input.  Two  distinct  classes  of  these  quantizers  exist;  those 
that  change  their  step  size  based  on  the  transmitted  value  of  the  error 
signal  and  those  that  change  their  step  size  based  on  the  variance  of 
the  unquantized  error  signal.  P.  Noll  * has  called  the  first  form 
"backward"  quantizers  since  they  look  backward  over  previously 
quantized  error  samples  to  adjust  their  step  size.  He  has  labelled 
the  second  form  "forward"  quantizers  because  they  look  forward  over 
the  unquantized  error  sample  to  obtain  their  step  size.  The  backward 
quantizers  need  not  send  the  quantizer  step  size  to  the  receiver 
because  the  receiver  can  regenerate  this  value  by  looking  at  the 
transmitted  sequence  representing  the  quantized  error  waveform.  This 
is  not  true  of  the  forward  quantizers.  Here,  since  the  step  size  is 
based  on  thevalue  of  the  unquantized  error  signal,  the  receiver 
cannot  regenerate  it  from  the  transmitted  sequence.  Thus,  forward 
quantizers  transmit  the  value  of  the  quantizer  to  the  receiver.  At 
a constant  data  rate,  coders  having  forward  quantizers  have  slightly 
lower  audio  input  bandwidth  than  do  those  having  backward  quantizers. 

Under  this  contract,  we  determined  that  both  forward  and 
backward  quantizers  can  be  made  to  work  equally  well  in  the  presence 
of  channel  errors.  This  contradicts  the  work  of  Noll  who  believed 
that  forward  quantizers  which  transmitted  their  step  size  would  work 
better.  At  the  time  of  his  writing,  however,  Noll  was  unaware  of 
modifications  to  the  backward  quantizers  which  would  improve  their 
tolerance  to  channel  errors . 

We  will  first  discuss  these  backward  quantizers  and  the 


modifications  needed  to  improve  performance  with  channel  errors. 
Then  we  will  describe  the  forward  quantizers  and  the  necessity  to 
send  the  step  size. 


19. 


2.2.1  The  Jayant  Quantizer 

The  Jayant  quantizer  adapts  its  step  size  instantaneously, 
expanding  or  contracting  its  va]ue  each  time  it  quantizes  a new  sample. 
Figure  2-5, illustrating  an  8 level  (3  bit/sample)  quantizer,  shows 
the  operation  of  this  quantizer.  If  the  incoming  sample  is  in  one 
of  the  inner  quantizing  levels,  the  quantizer  step  size  is  decreased 
for  the  next  sample.  If  the  incoming  sample  is  In  one  of  the  outer 
levels,  the  step  size  is  increased  instead.  These  changes  are  based 
on  the  quantized  output  and  thus,  the  Jayant  quantizer  is  classified 
as  backward  for  there  is  no  need  to  send  the  quantizer  step  size  as 
the  receiver  can  regenerate  it.  The  constantly  changing  step  size, 
which  is  attempting  to  track  the  input  signal  variance  is  given  by 

qk  = M(Ik)  qk-l 


at  time  k where  the  M(Ik)  represent  the  multiplication  factors  for 
the  step  size  on  the  kth  sampling  interval  and  1^  is  the  kth  transmissior 
symbol.  The  adaptive  quantizer  formulation  need  not  be  restricted  to 
quantizers  having  only  8 levels.  Any  quantizer  having  more  than  2 levels 
can  be  made  adaptive  with  the  Jayant  algorithm.  In  fact,  a 3 level 
quantizer  Is  possible.  This  quant izer , which  uses  (1.6  bits/sample) 
quantizer,  expands  the  step  size  whenever  its  two  outer  values  are  used 
and  decreases  the  step  whenever  the  inner  value  is  used.  The  three 
level  quantizer  is  desirable  whenever  the  data  rates  must  not  exceed 
9600  bps,  whereas  the  8 level  quantizer  requires  at  least  16000  bps  at 
2500  Hz  bandwidth. 

1 
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In  general,  for  a given  transmission  rate,  there  is  an 
inverse  relationship  between  the  number  of  quantizer  levels  and  the 
audio  input  bandwidth.  At  16,000  bps  an  8 level  quantizer  (3  bits/ 
sample)  limits  the  audio  bandwidth  to  about  2500  Hz,  while  a 5 level 
quantizer  (2.33  bits/sample)  would  permit  about  3200  Hz  bandwidth 
and  a ^ level  quantizer  would  allow  more  than  3300  Hz.  Thus,  using 
a 5 or  6 level  quantizer  instead  of  an  8 level  quantizer,  the  audio 
bandwidth  of  the  coder  will  increase  from  2500  Hz  but  the  quantizing 
noise  will  also  increase  because  fewer  quantizing  levels  are 
employed . 

The  Jayant  algorithm  performs  well  without  channel  errors. 
Unfortunately,  channel  errors  significantly  impair  performance.  To 
see  this  consider 


q 


k 


‘Sk-l 


k 

11  M(Ii)  % 

1 = 1 1 


where  1^  is  the  ith  transmitted  character  and  MCI^)  is  the  multiplier 
corresponding  to  1^ . 

Assume  now  that  one  transmission  error  is  made  at  time  1 

causing 

M(I^)  to  become  M(I^') 


qz  to  become  q^'  = MCI^')  qJl_1 


and 
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afterwards  the  new  quantizer  levels  for  time  greater  than  £.  becomes 

q ' =■  S N(I,)  Mtr,')  qo  • 

L=1 
Lfil 

or  simplifying  we  obtain 

, m(i4') 

qk  = qk  M(I4) 

Consequently,  the  effects  of  even  a single  channel  error  never  die 
out.  In  fact,  an  error  causes  a level  shift  in  qk  by  the  factor 
M(I4' )/M(I4) . Simulations  and  real-time  implementation  show  that 
errors  cause  noticeable  fades  and  increases  in  output  level.  Only 
when  the  quantizer  saturates  at  its  maximum  level  or  cuts  off  at 
its  minimum  level  does  the  quantizer  again  track  properly.  Thus, 
the  Jayant  algorithm  must  be  altered  if  it  is  to  be  useable  over 
real  communication  channels. 


M(IJ 
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2.2.2  The  Modified  Jayant  Algorithm 


Wilkinson  and  Goodman11  have  recently  shown  how  a simple 
change  to  the  Jayant  algorithm  can  make  this  algorithm  robust  to 
channel  errors.  Their  modification  was 

k k-i 

qk  = M(Ik)qk-i  a = 1;1  "(V®  % 

where  1^  is  the  ith  transmitted  symbol  and  M(I^)  is  the  multiplier 
corresponding  to  1^.  Again  if  one  channel  error  curves  at  time  £ 

M(I£)  - MCI^  ) 


and 


k 

1 

<-i  , k-£ 

M( If  ) 

qi/  = * 

i=l 

Md^l0 

”<*£  > qQ  * 

i = £ 

simplifying 

r~  • ■"n 

k-£ 

1 

) 

a 

qk  = qk  ' 

M(I.) 

k-£ 


Assuming  a < 1,  then  for  times  k far  removed  from  £ (k  >>  £) 


qk  * qk 


The  speed  of  this  convergence  depends  on  how  near  a is  to  1.  If  £ 
is  not  too  close  to  unity,  convergence  is  rapid. 

Consequently,  the  effects  of  a single  channel  error  decay 
with  time  for  the  modified  Jayant  algorithm  and  this  algorithm  is 
much  more  suitable  for  use  over  communications  channels. 
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2.2.3  The  Forney  Quantizer 

The  Forney  quantizer  was  another  attempt  at  creating  a 
backward  quantizer  which  is  Insensitive  to  channel  errors. 

As  shown  in  Figure  2-6, 


_ - 
qk  ' 2 


ek  = Gk  + Ck 


where 


and 

°k  ls  “Vi + Hi(Ik> 

Ck  = eck-l  + M2^k^ 
where  specific 

M^( 1^)  M2(Ik)  a & 6 are  given  in  Figure  2-6. 

Consequently,  this  quantizer  also  increases  its  step  size  when  large 
inputs  occur  and  contrasts  the  step  size  during  low  level  Inputs. 

The  quantizer  step  size  is  non-uniform  and  reacts  rapidly  to 
occasional  large  spikes  caused  by  pitch  pulses. 

Noting  that 

Gk  = 2 ^ 1 Ml(Ii} 

K i=l  11 


and 


Ck  = Z 
K i=l 


„k-i 


m2(i1) 


we  have 
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Figure  2-6.  Forney  Adaptive  Quantizer 


at  time  £ assume  a channel  error  causes 


then 


W * «2(It  > 


2/^1  ak'1M1(I1)  + ak-^M^(  ' ) + J6k_iM2(I1)  + 6k"lM2<It' 


by  adding  and  subtracting 


ak_*Mi(I£)  and  ek‘2,M2(i]l) 


to  the  exponent 
we  have 


as  k >>  £ -*•  0 

for  a < 1 
and  qk'  - qk 

Thus  the  Forney  quantizer  also  is  insensitive  to  channel  errors 
because  after  a channel  error  the  step  size  approaches  the  desired 
step  size  with  time.  Moreover,  this  algorithm  responds  more  rapidly 
to  sudden  changes  in  the  input  than  does  the  modified  Jayant 
algorithm. 


r 
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2.2.4  The  Fixed/Frame  Quantizer 

1 2 

Fixed/frame  quantizers  as  originally  used  in  APC  algorithms  ’ 

calculate  a step  size  for  an  entire  frame  of  data.  Since  the  receiver 

cannot  regenerate  this  data  from  the  quantized  error  signal  sequence, 

the  step  size  must  be  sent  to  the  receiver  once  each  frame.  This  is  a 

q in 

forward  type  of  quantizer  as  described  by  Noll.  ’ 

The  step  size  Q is  a function  of  the  RMS  energy  in  the  input 
signal.  The  energy  U is  calculated  during  estimation  of  the  predictor 
coefficients  and 


Q = 1/2 (U/T )1/2 


where  T = the  frame  length  and  the  factor  1/2  was  determined 
experimentally. 

The  fixed/frame  quantizer  is  insensitive  to  errors  because  a 
single  channel  error  effects  only  one  value  of  the  reconstructed  error 
signal.  Of  course,  if  an  error  occurs  in  transmissing  Q then  an  entire 
frame  of  data  is  shifted  in  level.  Fortunately,  the  following  frames 
of  data  are  unaffected. 


;*y . 

2.3  Fortran  Fixed- Point  Simulations 


The  various  APCQ  configurations  and  error  signal  quantizers  were 
investigated  through  FORTRAN  fixed-point  simulations  intended  to  implement 
the  algorithms  with  the  same  numerical  accuracy  as  the  real-time  software. 
Writing  these  codes  In  FORTRAN,  perm  1 tted  us  to  Investigate  many  designs 


rapidly  without  the  long  debugging  and  check  out  times  associated  with 
writing  real-time  software. 

The  block  diagram  of  the  APCQ  coder  is  shown  in  Figure  2-7 
with  a more  detailed  indication  of  the  processing  at  the  analyzer 
(transmitter)  shown  in  Figure  2-3-  Figure  2-8  indicates  the  general 
flow  of  data  in  the  analyzer  with  Figure  2-9  indicating  the  data 
flow  of  the  synthesizer  (receiver).  These  flowcharts  indicate  that  the 
speech  stored  on  digital  tape  is  read  in  a frame  at  a time,  processed 
by  the  fixed-point  algorithms,  and  reconstructed  and  placed  on  an  output 
tape  for  listening.  The  I/O  is  handled  by  the  digital  tape  handlers 
written  by  GTE  Sylvania  prior  to  this  contract.  Thus  the  I/O  for  the 
simulation  was  straightforward,  and  effect  was  placed  on  the  scaling 
and  normalization  needed  to  perform  the  processing  operations  in  fixed- 
point  arithmetic.  Appendix  A provides  a complete  listing  of  the  fixed- 
point  simulation  software  with  flow  charts  indicating  the  scaling  and 
normalization  operations  in  each  arithmetic  function. 

2.^  Simulation  Results 

2.^.1  Coder  Performance  without  Channel  Errors 

All  of  the  speech  coders  of  this  contract  attempted  to  preserve 
the  time  waveform  shape.  As  a result,  a reasonable  performance  measure 
for  these  systems  is  signal-to-noise  ratio  (S/N)  as  defined  in  Figure 
2-10  Informal  listening  tests  of  these  coders  indicate  that 
listeners  generally  prefer  those  systems  having  higher  S/N  ratios  over 
those  having  lower  ratios.  The  computer  simulations  developed  S/N 
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8-LEVEL  QUANTIZER 


Figure  2-7,  Simulation  of  APCQ  Speech  Encoding  System 
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Figure  2-8.  Analyzer  Flowchart  of  Fixed-Point  Simulation  for  APCQ  Coder 
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ratios  for  each  non-overlapping  25  msec  speech  segment.  These  ratios 
were  then  averaged  to  obtain  a cumulative  S/N  ratio  for  an  entire 
sentence.  It  is  these  cumulative  S/N  ratios  that  are  used  In  our 
results  presented  in  the  graphs  of  Figures  2-11  through  2-12. 

Figure  2-11  shows  the  performance  of  APCQ  systems  using  an  8 level 
Jayant  quantizer  operating  at  16  kbps.  This  data  was  obtained  by 
averaging  the  cumulative  S/N  ratio  obtained  for  each  sentence  over  a 
total  of  six  input  test  sentences  containing  both  male  and  female 
voices.  Each  test  sentence  contained  about  100  analysis  frames  and 
thus,  each  point  on  Figure  2-11  represents  the  average  of  600  frames. 
Referring  to  Figure  2-11,  the  performance  increases  with  increasing 
predictor  order  and  with  the  inclusion  of  a pitch  loop.  The  zeroth 
order  predictor  without  pitch  is  equivalent  to  APCM.  This  system 
performs  roughly  2-3dB  worse  than  the  CVSD  system  we  have  chosen  as 
our  benchmark.  Other  simulations  performed  as  part  of  this  study 
indicate  that  ADPCM  has  a S/N  equivalent  to  APCQ-1  (APCQ  with  1 
predictive  coefficient)  without  pitch.  As  Figure  2-11  indicates, 

APCQ-1  (without  pitch)  is  roughly  equivalent  to  CVSD.  The  addition 
of  a pitch  loop  adds  between  1.5  and  2.0  dB  to  the  S/N  performance. 

2.4,2  Coder  Performance  with  Channel  Errors 

Figure  2-12  shows  the  performance  of  these  coders  in  the  presence 
of  channel  errors  for  the  Jayant,  Modified  Jayant,  Forney  and  fixed 
frame  quantizers.  Controlled  random  bit  errors  were  applied  to  all 
transmitted  data  equally.  The  coder  used  to  evaluate  these 
quantizers  was  a first  order  adaptive  predictor  (APCQ-1)  without 
a pitch  loop.  In  preparing  the  data  of  this  Figure,  only  one  test 
sentence  was  used  due  to  large  amount  of  processing  involved.  Using 
selected  points  taken  from  other  test  sentences  showed  similar 
trends  in  performance. 


Figure  2-11:  S/N  of  CVSD  and  Adaptive  Predictive  Coders  with 

Jayant  Quantizer  [APCQ]  vs.  Predictor  Complexity 


COEFFICIENT  PREDICTOR  ORDER 


As  can  be  seen  from  Figure  2-12,  the  Jayant  quantizer,  CVSD  and 
the  Forney  algorithms  are  within  ldB  of  each  other  in  the  absence  of 
channel  errors.  Under  this  condition,  the  fixed/frame  quantizer  is 
3dB  below  the  performance  of  these  other  quantizers  and  the 
modified  Jayant  quantizer  is  ldB  below  the  fixed/frame  quantizer. 

The  Jayant  algorithm,  however,  is  extremely  sensitive  to  channel 
errors.  Its  performance  falls  off  rapidly  when  errors  occur.  This 
quantizer  produces  noticeable  distortion  at  error  rates  as  low  as 
10  and  is  unacceptable  at  BER's  of  10  J and  above.  The  Forney 

quantizer  is  more  robust,  maintaining  good  performance  at  10  , 

-2  —1 
fair  performance  at  10  and  unacceptable  performance  at  10  . While 

starting  at  a lower  performance  level  than  the  other  quantizers,  the 

fixed  frame  quantizer  does  not  degrade  rapidly  with  increasing 

channel  errors.  In  fact,  it  provides  better  performance  than  either 

_2 

the  Jayant  or  Forney  quantizers  at  error  rates  higher  than  10 
The  modified  Jayant  quantizer  using  a decay  factor  of  O.96  had 

a S/N  that  was  always  below  the  fixed  frame  quantizer,  but  at  error 

_2 

rates  above  10  , was  higher  than  the  FORNEY  quantizer.  Note  that 

CVSD,  the  simplest  system,  exhibited  the  best  performance  against 
channel  errors  until  error  rates  exceed  10_\ 

Based  on  the  simulation  data  of  Figure  2-12,  the  fixed  frame 
quantizer  was  chosen  to  evaluate  the  performance  of  the  APCQ  system 
with  more  complex  predictors.  Figure  2-13  shows  how  APCQ  systems 
using  a fixed  frame  quantizer  perform  against  channel  errors  over  the 
same  sentence  material  as  Figure  2-12.  In  this  figure,  the  channel 
errors  are  imposed  uniformly  on  the  transmitted  data  including  the 
pitch  M,  pitch  gain  a,  the  PARCOR  coefficients  and  the  quantizer 
step  size  as  well  as  on  the  quantized  error  signal  sequence. 
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At  low  error  rates,  a fourth  order  APCQ  system  with  pitch 
outperforms  CVSD  by  3dB.  As  bit  errors  are  Imposed,  the  APCQ 
system  loses  its  advantage  over  CVSD.  At  high  error  rates,  the 
presence  of  the  pitch  loop  becomes  detrimental  to  APCQ  performance. 

At  error  rates  between  10-^  and  10~^,  the  APt)Q  system  with  fixed/frame 
quantizer  and  CVSD  have  similar  S/N  rates.  Note  that  the  simplest 
system  (i.e.,  APCM)  outperforms  CVSD  at  error  rates  of  10  ^ . This  is 
shown  by  the  bottom  curves  of  Figure  2-13  with  predictor  order  equal 
to  zero  and  no  pitch  loop  present.  Figure  2-14  provides  an  alternate 
method  of  graphing  some  of  the  data  given  in  Figure  2-13  • These 
curves  compare  the  channel  error  performance  of  a fourth  order  adaptive 
predictor,  with  and  without  pitch,  to  CVSD.  At  low  error  rates, 
fourth  order  systems  are  superior  to  CVSD.  At  higher  error  rates, 

CVSD  is  either  superior  to  APCQ  with  pitch  or  comparable  to  APCQ 
without  pitch. 

These  results  indicate  that  system  performance  is  a critical 
function  of  channel  error  rate.  As  low  errors  occur,  the  complex 
APCQ  systems  outperform  CVSD  at  16  Kb/sec.  For  high  channel  error 
environments,  CVSD  is  superior.  If  high  quality  is  the  objective  in 
low  error  environments,  then  the  APCQ  system  provides  the  designer 
with  a superior  voice  digitization  technique. 

Based  on  these  simulation  results,  GTE  Sylvania  and  DCA 
agreed  to  develop  a final  system  using  a fixed  frame  quantizer.  To 
increase  the  audio  bandwidth  from  2500  Hz  to  3200  Hz,  the  8 level 
quantizer  was  reduced  to  a 5 level  quantizer.  This  lowered  the  S/N 
ratio  of  the  APCQ  system  approximately  3dB  from  those  presented 
in  Figures  2-11  through  2-l4. 

2.4.3  Improved  Error  Signal  Quantizers 

The  modified  Jayant  algorithm,  however,  did  not  become 
available  until  after  we  made  this  decision  to  implement  with  the 
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Figure  2-14.  S/N  of  CVSD  and  4th  Order  Adaptive  Predictive  Coder  (Fixed 
Quantizer)  vs.  Bit  Error  Rate  @ 16  KBPS 
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fixed/frame  quantizer.  Moreover  as  Figure  2-12  illustrates,  its  S/N 
is  always  lower1  than  the  other  quantizers.  Nevertheless,  Informal 
listening  tests  conducted  after  the  decision  to  implement  with  the 
fixed/frame  quantizer  indicated  little  if  any  degradation  between 
the  Jayant  and  modified  Jayant  quantizers  without  channel  errors. 

Even  more  important  was  the  discovery  that  the  distortions 
introduced  by  the  modified  Jayant  quantizer  were  less  objectionable 
than  those  introduced  by  the  fixed/frame  quantizer.  Consequently, 
during  the  last  weeks  of  the  contract,  additional  simulation  runs 
were  performed  to  observe  the  behavior  of  the  modified  Jayant  algorithm 
further . 

Figure  2-15&16  illustrates  the  performance  of  the  APCQ  coder 
with  channel  errors  inserted  into  the  transmitted  data  using  the 
modified  Jayant  quantizer.  The  S/N  ratios  behave  about  the  same  way 
as  the  fixed  frame  quantizer  in  the  presence  of  channel  errors.  Thus 
for  error  exponents  greater  than  10  , the  systems  having  no  pitch 

predictor  and  a low  order  coefficient  predictor  performed  better  than 
the  systems  with  a pitch  loop  and  high  order  predictions.  The  simple 
systems  are  better  since  the  channel  errors  occur  on  the  pitch 
parameters  and  predictor  coefficients  and  these  must  significantly 
affect  the  S/N  ratio.  On  good  channels  (BER  < 10~^)  the  S/N  ratio 
improves  as  the  predictor  becomes  more  sophisticated. 

The  S/N  ratio  of  the  APCQ  coder  with  modified  quantizer  is 
always  lower  than  that  for  the  fixed/frame  quantizer.  The  perceived 
voice  quality,  however,  is  generally  higher  because  the  types  of 
distortion  are  different. 

Thus  the  S/N  curves  while  indicative  of  the  expected  performance 
of  a technique  are  not  always  directly  related  to  voice  quality. 


20 


S/N  OF  ADAPTIVE  PREDICTIVE  CODERS  USING  8 LEVEL 
MODIFIED  JAYANT  QUANTIZER 


Figure  2-15.  S/N  of  Adaptive  Predictive  Coders  Using  8-Level 
Modified  Jayant  Quantizer 
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BIT  ERROR  RATE 


Figure  2-16.  S/M  of  CVSD  and  4th  Order  Adaptive  Predictive  Coder 

(Modified  Jayant  Quantizer)  versus  Bit  Error  Rate  at  16  Kb/s 
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The  modified  Jayant  quantizer  has  high  voice  quality  but  fairly 
low  S/N  ratios. 

The  modified  .Jayant  quantizer  because  It  does  not  degrade  as 
the  predictor  becomes  less  complex  offers  the  potential  of 
simplifying  the  APCQ  system  dramatically. 
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SECTION  3 

Real-Time  Software 
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3.1  Description  of  APCQ  Program 

3.1.1  General  Discussion 

The  real-time  APCQ  program  written  1 ri  assembly  language  enables 
the  EDM  to  collect,  analyse,  transmit,  receive  and  syntheize  speech  data 
in  full  duplex  mode.  The  sequence  of  operations  involves  initialization 
of  the  system  parameters,  set-up  of  the  interrupt  processing,  acquisi- 
tion of  synchronization,  processing  of  speech  via  the  APCQ  algorithm, 
transmission  and  reception  of  parameters  and  reconstruction  of  speech. 

The  APCQ  algorithm  implemented  is  a fourth  order  adaptive  pre- 
dictive coder  with  a pitch  loop  and  a five- level  fixed  frame  quantizer. 
Each  speech  sample  is  predicted  as  a weighted  linear  combination  of  four 
past  speech  samples.  The  error  signal,  computed  from  the  difference  of 
the  actual  and  the  estimated  speech  sample,  is  then  quantized  into  five 
levels  by  the  quantizer  shown  in  Figure  3-1  with  level  estimated  from 
the  RMS  power  of  the  error  signal  which  remains  constant  throughout  a 
whole  frame  of  speech  data. 

The  key  to  correctly  sequencing  the  operations  is  the  use  of  two 
interrupts  in  the  EDM:  the  speech  side  interrupt  and  the  line  side 

interrupt.  The  speech  side  interrupt  is  set  to  interrupt  the  computer 
at  regular  intervals  and  transfer  control  to  software  (starting  at 
Location  0),  which  inputs  one  speech  sample  into  the  computer  and  out- 
puts one  speech  sample  from  the  computer.  Careful  count  is  kept  of  the 
number  of  samples  read  in;  and  when  a fra m'-  of  data  haa  b<-<-ri  Input,  an 
indicator  (transmit  flag)  is  set  showing  the  analysis  can  begin. 

The  line  side  interrupt  also  regularly  Interrupts  the  computer 


and  transfers  control  to  software  (starting  at  Location  1),  which  out- 
puts modem  clock  (five-volt  levels  from  a flip-flop)  and  one  bit  of 
transmission  data.  At  every  eighth  interrupt  it  inputs  eight  bits  of 
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received  data  which  have  been  accumulating  in  an  eight-bit  shift 
register.  When  it  receives  a frame  of  input  data,  it  sets  an  indicator 
(receiver  ready  flag)  showing  that  synthesis  can  begin. 

Since  the  program  operates  in  less  than  real-time,  it  executes 
instructions  in  a wait  loop  while  it  idles.  It  leaves  the  wait  loop  to 
analyze  speech  whenever  it  sees  that  the  transmit  ready  flag  has  been 
set  by  the  speech  interrupt.  After  analysis  it  returns  to  the  wait  loop 
where  it  will  leave  again  for  the  synthesis  when  the  line  side  inter- 
rupt raises  the  receive  ready  flag.  In  addition,  it  also  interrupts  the 
sync  flag  and  leaves  for  the  sync  search  routine  whenever  the  flag  is 
set.  Naturally,  the  speech  and  line  interrupts  are  functioning  even 
during  execution  of  the  wait  loop.  Speech  is  input  and  output  at  each 
sampling  instant,  and  data  is  sent  and  received  at  16,000  bps. 

3.1.2  Initialization  and  Wait  Loop 

When  the  program  is  started,  it  immediately  goes  to  an  initial- 
ization mode  where  all  parameters  and  counters  of  analyser  and  synthe- 
sizer filters  are  cleared.  It  then  goes  into  the  wait  loop  shown  in 
flow  chart  Figure  3-2  and  waits  until  one  of  three  flags,  namely: 
transmitter  ready  flag,  receiver  ready  flag  and  sync  search  ready  flag 
is  set  where  it  performs  the  corresponding  processing.  At  the  end  of 
the  operation,  the  program  returns  to  the  wait  loop  to  use  up  the  idle 
time  left  over.  Reinitialization  can  always  be  achieved  by  setting  one 
of  the  front  panel  switches. 

3.1.3  Synchronization 

After  received  80  bauds  (1.5  second)  of  data,  the  sync  search 


I 

I 


ready  flag  is  set  where  the  program  goes  to  sync  search  loop  shown  in 
Figure  3-3. 
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Figure  3-3.  Synchronization  Search 


To  understand  the  synchronization  process,  it  is  necessary  to 
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describe  the  format  of  the  transmission  data.  Figure  3-4  shows 
that  the  transmission  baud  consists  of  312  bits  arranged  in  the 
following  manner:  280  bits  for  the  predictor  error  signal,  one 

synchronization  bit  which  is  always  set  to  a one,  three  bits  for  fourth 
PARCOR  coefficient,  three  bits  for  the  third  PARCOR  coefficient,  four 
bits  for  the  second  PARCOR  coefficient,  four  bits  for  the  first  PARCOR 
coefficient,  three  bits  for  the  pitch  gain  ALPHA,  six  bits  for  the  pitch, 
four  bits  for  the  RMS  power  of  error  signal,  three  bits  for  the  normalizer 
and  a second  one  bit  for  synchronization  which  is  also  always  a one. 

The  receiver  develops  a histogram  of  the  received  data,  counting  the 
times  a one  is  received  in  each  of  the  312  bit  positions.  After  80 
frames  (approximately  1.5  seconds)  the  receiver  examines  the  histogram 
for  positions  of  two  values  near  80  (the  values  may  be  slightly  less 
than  80  because  of  transmission  errors)  which  are  31  bit  positions 
apart.  From  the  position  of  the  sync  points,  the  software  then  knows 
where  the  last  bit  in  the  frame  is  located,  and  it  shifts  the  received 
sequence  so  that  the  first  bit  of  the  baud  lines  up  with  the  first 
storage  location.  Once  synchronization  has  been  achieved,  the  locations 
of  the  synchronization  bits  are  saved  and  compared  with  that  of 
subsequent  sync  searches  to  determine  the  necessity  of  shifts  in  the 
following  received  sequence. 

3.1.4  Line-side  Interrupt 

Figure  3-5  describes  the  program  flow  whenever  the  clock 
interrupts  the  computer  for  data  transmission  or  reception  by  trapping 
to  Location  1.  First  the  flip-flop  which  drives  the  modem  clock  is  set 
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Figure  3-4.  Format  of  Transmission  Data 
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Figure  3-5.  Line  Side  Interrupt 
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or  reset  depending  on  whether  the  interrupt  was  even-numbered  or 
odd-numbered.  Whenever  it  is  set,  one  bit  of  data  is  also  sent  and 
the  receiver  buffer  is  interrogated  to  see  if  eight  bits  of  data  have 
been  received.  If  received  data  is  ready,  the  software  inputs  the 
data,  compiles  the  histogram,  and,  whenever  the  baud  is  complete, 
sets  a receiver  ready  flag.  Every  80  bauds  it  sets  a flag  indicating 
that  a synchronization  search  could  begin  if  desired. 

8.1.5  Speech  Side  Interrupt 

The  speech  side  interrupt  is  straightforward , as  shown  in 
Figure  3-6.  When  the  computer  interrupts  and  traps  to  Location  0, 
it  reads  in  one  speech  sample,  reads  out  one  sample,  and  when  the 
buffer  is  full,  sets  the  transmit  ready  flag.  The  speech  data  is  double- 
buffered  to  prevent  loss  of  data  during  analysis. 

3.1.6  Speech  Analysis  and  Synthesis 

Flow  charts  of  Figures  3-7,  3-8,  3-9  and  3-10  describe  the 
operations  of  the  analysis  of  the  APCQ  coder.  Upon  entry  to  the  analyser 
the  input  data  is  first  normalized.  Then  pitch  is  computed  via  the  AMDF 
function  as  shown  in  Figure  3-5.  Next  the  pitch  gain  ALPHA  and  the  reduced 
waveform  are  calculated  as  illustrated  in  Figure  3-7. 

In  Figure  3-8  the  autocorrelation  coefficients  of  the  reduced 
waveform  are  determined  that  are  later  used  to  compute  PARCOR  coef- 
ficients via  the  matrix  inversion  routine  shown  in  Figure  3-9  from 
which  actual  predictor  coefficients  (a's)  are  derived  as  in  Figure 
3-14.  Finally,  the  RMS  power  of  the  error  signal  is  calculated  and 
used  to  quantize  error  signal  and  coded  Into  seven  bit  words  as  shown 
in  Figure  3-10. 

Figure  3-12  and  3-13  describes  the  operations  of  the  synthesizer 
where  all  transmitted  parameters  are  decoded  and  speech  is  reconstructed 
by  inverse  filtering. 
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Figure  3-6.  Interrupt  Loop  Speech  Side 
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Figure  3-7.  Analyzer  for  APCQ  Coder 
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Figure  3-11.  Transmission  Data  Construction 
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Figure  3-12.  Data  Reception,  Decoding  and  Dequantization 
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Figure  3-15.  AMDF  Pitch  Calculation 
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3.2  Operation  Procedures  for  the  APCQ  Program 

3.2.1  Connecting  the  System 

a.  Connect  the  EDM's  to  each  other  either  in  back-to-back 
mode  shown  In  Figure  3-17  or  through  a 16,000  bps  modem 
as  shown  in  Figure  3-]  8. 

1.  Use  the  terminal  marked  "Date.''  on  the  rear  of  the  CPU  to 
connect  to  the  modem  or  the  other  EDM. 

2.  Place  the  toggle  switch  to  appropriate  RS-232  or  MIL-188 
interface  position  as  required 

3.  Strap  the  modem  to  operate  from  the  external 
transmit  clock  or  the  leading  edges  of  the  data 

b.  Connect  the  handsets  or  audio  tape  recorders  to  the 
proper  input  and  output  jacks  at  rear  of  CPU 

3.2.2  Read!  ig  in  the  EDM  Software 

a.  Place  all  switch  register  toggle  switches  in  the 
up  position 

b.  Power  on  the  PSP  and  card  reader 

c.  Depress  the  "PM  CLEAR"  switch  to  clear  program  memory 

d.  Set  "NORMAL/TEST  PMIR"  switch  to  "NORMAL" 

e.  Set  "RUN/STOP"  switch  to  "STOP” 

f.  Load  the  binary  deck  and  depress  the  "MOTOR"  and 
then  "START"  switches  on  the  card  reader 

g.  Cards  (load  PM,  load  DM,  PMMA  checksum,  DM  checksum) 
are  read  in.  Card  reader  halts  before  it  gets  to  PMIR 
checksum  cards 

h.  Depress  "STOP"  switch  on  card  reader 


'igure  3-17.  System  Test  Configuration 


Modem  Interface  Signals  and  Interconnection 


I.  Set  "NORMAL/TEST  PMIR"  switch  to  "TEST  PMIR . " 
Depress  "MANUAL  INST"  switch 

J . Restart  the  card  reader  by  depressing  "MOTOR"  and 
"START"  switches 


k.  After  the  whole  deck  is  read,  set  "NORMAL/TEST  PMIR" 
switch  back  to  "NORMAL" 

l.  Set  "INTERRUPT  PROGRAM/DISABLE"  switch  to  "INTERRUPT 
PROGRAM" 


m.  Set  "RUN/STOP"  switch  to  "RUN"  to  start  executing 
the  program. 


3.?.  3 Initializing  the  EDM 

After  the  program  has  been  running,  reinitialization 
can  be  achieved  by  simultaneously  holding  switch  register 
toggle  switch  15  In  the  down  position  for  each  EDM  and  then 
returning  it  to  the  up  position. 
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3.3  Description  of  Continuous  Variable  Slope  Delta  Modulation  Program 
3.3.1  Basic  Principles  of  Operation 

CVSD  represents  one  of  the  simplest  voice  digitizers  known, 
since  there  are  relatively  few  mathematical  operations  per  input  sample. 
Figure  3-19  depicts  a digital  implementation  of  a CVSD  transmitter  and 
receiver  operating  on  the  DCA  EDM's.  The  input  speech  is  first  con- 
verted to  digital  by  a 12  bit  A/D  converter  after  lowpass  filtering 
through  an  elliptic  filter  having  a 3 dB  frequency  of  3200  Hz.  The  CVSD 
algorithm  subtracts  an  estimate  of  the  incoming  sample  from  the  input 
and  quantizes  the  difference  or  error  sample  to  +1  if  the  error  is  posi- 
tive and  -1  if  negative.  The  quantized  value  is  encoded  for  transmis- 
sion by  the  mapping  1 into  1 and  -1  into  0.  The  algorithm  then  uses  a 
two  state  process  to  construct  the  estimate  for  the  next  Incoming  sample. 
The  first  stage  computes  a quantizer  gain  factor  by  interacting  with 
time  constant  T^  and  uses  this  factor  to  adjust  the  level  of  the 
quantized  error  signal.  The  second  stage  filters  this  new  error  signal 
through  an  integrator  having  time  constant  T2  to  produce  the  new  speech 
estimate.  The  CVSD  receiver  then  reconstructs  the  +1  bit  stream  of  the 
transmitter's  error  quantizer.  This  data  excites  the  same  quantizer 
gain  and  error  signal  integrator  (which  formed  the  speech  estimate  at 
the  transmitter)  and  generates  the  synthesized  speech  after  D/A  conver- 
sion and  low  pass  filtering. 

With  varying  gain  in  the  error  signal,  the  system  adapts  to 
the  slope  of  incoming  speech  waveform.  Severe  slope  overload  occurs 
with  insufficient  quantizer  gain  perceived  as  muffled  speech  with  drop 
of  volume.  On  the  other  hand,  excessive  gain  will  cause  overshoots 
which  introduce  granular  noise  into  the  synthesized  speech.  In  this 
study,  three  slope  adaptation  stretegies  are  investigated: 
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Figure  3-19.  Transmitter  and  Receiver  of  CVSD-B 
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3 • 3 • ^Standard  Three  Bit  Memory  (CVSD-B) 

Figure  3-19  shows  the  block  diagram  of  the  CVSD-B  transmitter 
and  receiver.  By  employing  a three  bit  shift  register,  memories  of 
the  present  and  past  two  error  samples  are  retained  and  utilized 
to  estimate  the  incoming  speech.  Whenever  the  three  bits  of  the 
shift  register  are  identical  (111  or  000),  a waveform  with  steep 
slope  is  usually  present.  Then  a spike  of  amplitude  A is  introduced 
into  a single  pole  smoothing  filter  with  coefficient  a^.  Whenever 
the  three  bits  are  different,  the  logic  introduced  a zero  into  the 


same  filter  shown  as: 


f(J)  = a,  f(j-l)  + 


A , if  shift  register  = 000  or  111 
0 , otherwise  (3-1) 


Then  the  filter  output  for  a constant  input  amplitude  A or  0 is 


given  by 


f(n)  = 


(1-a,  1 ) ; if  shift  register  * 000  or  111 


; otherwise 


(3-2) 


When  f(j)  is  added  to  a small  bias  6,  the  result  is  a gain  factor 
for  multiplying  the  +1  data  coming  from  the  error  quantizer.  This 
new  waveform  then  generates  the  speech  waveform  by  exciting  a second 


single  pole  integration  filter  with  coefficient  a^. 

In  order  to  determine  the  amplitude  of  A that  yields  the  best 


system  performance,  a parameter  called  compression  ratio  (CR)  which  is 
defined  as 


CR  = f max  + 6 
6 


where  f 


max  1-a, 


(3-3) 


was  investigated.  This  compression  ratio  is  related  to  A by 


<5*  (CR-1 ) 


(3-H) 


i 


compression  ratio,  the  system  tracks  the  incoming  waveform  with  a 
smaller  slope  overload. 


As  shown  In  Eq.  ( 3-1*  ),  high  pulse  height  A can  also  be  obtained 
by  raising  the  bias  value  6.  Hence  a high  5 tends  to  yield  a better 

estimate  to  a fast  changing  input  waveform.  But  in  the  presence  of 
a slow  varying  waveform,  such  as  silence  intervals  between  speech, 
the  main  excitation  component  to  the  integration  filter  is  6 which 
yields  an  oscillatory  reconstructed  waveform  with  a value  +6. 

Hence  large  6 tends  to  increase  the  granularity  of  the  processed 
speech . 

Besides  the  above  parameters,  the  time  constants  and 
of  the  gain  and  prediction  filters  also  contribute  to  the 
reconstruction  of  speech  waveform.  A large  value  for  the  time 
constant  T^,  which  is  related  to  the  coefficient  a1  by 

Tx  * -T/ln  ^ (3-5) 

Where  ^ is  the  transmission  rate  tends  to  have  a smoothing  effect 
on  producing  the  quantizer  gain  which  decreases  the  sensitivity 
to  channel  errors.  A small  value  for  time  constant  which  is 
given  by 

T2  = -T/ln(a2)  (3-6) 

produces  a more  energetic  tracking  of  the  input  waveform. 


I 

I 


exponentials  is  shown  in  Figure  3-20.  This  algorithm  uses  a 3-bit 
shift  register  with  the  last  bit  corresponding  to  the  present  error 


sample.  The  content  of  this  shift  register  are  used  to  generate  pulses 
of  different  length  which  excite  a smoothing  filter.  Table  3-1  provides 
the  amplitudes  of  these  pulses  as  a function  of  the  shift  register 


contents . 

Content  of  Shift  Register-  Value 

000  3/16 

001  0 

010  -5/128 

011  5/64 

100  5/64 

101  -5/128 

110  0 

111  3/16 


Table  3-1:  Pulse  Heights  for  FORNEY  CVSD  Quantizer 

The  output  of  this  filter  after  addition  of  a small 

x 

bias  6,  is  used  as  the  exponent  on  the  expansion  2 . For 
implementation  purposes,  an  approximation  to  the  exponential 
expression  is  given  by 

2(IQK  ♦ «KF),  2(IQK)  (1+qkf)  (3_7) 

where  IQK  = integer  part  of  the  filter  output  and  bias 

QKF  = fractional  part  of  the  filter  output  and  bias 
The  result  of  this  exponentiation  is  then  used  to  excite  an 
integration  filter.  Because  exponentials  are  used  for  both 
increasing  and  decreasing  the  quantizer,  a faster  gain  can  be 


1 


1 


[ 


I 

I 

I 
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achieved  within  a relatively  short  time,  resulting  in  a better  dynamic 
range  of  the  system  and  better  tracking  of  rapidly  increasing  and 
decreasing  waveforms. 

3.3-4  Jay ant's  One  Bit  Memory  (CVSD-D) 

Figure  3-21  illustrates  the  Jayant  slope  adaptation  logic. 

In  contrast  to  the  two  previous  algorithms,  the  Jayant  method 
uses  only  one  bit  of  memory.  Every  time  an  error  sample  is  generated, 
it  is  compared  with  the  previous  quantized  value.  If  they  are  differ- 
ent, the  present  quantizer  level  is  multiplied  by  a factor  Q < 1.  If 
they  are  the  same,  it  implies  a continuous  rising  or  falling  waveform. 
Then  the  present  level  is  multiplied  by  a factor  P > 1.  They  are  used 
to  excite  the  integration  filter  to  obtain  an  estimate  of  the  incoming 
speech  waveform. 

3.4  Discussion  of  Real-Time  Code 

The  real  time  software  of  CVSD  is  straightforward  when  com- 
pared to  APCQ  program  codes.  The  sequence  of  operations  involve  the 
initialization  of  speech  and  line  side  interrupts,  speech  processing 
using  one  of  the  quantizer  gain  update  strategies  as  discussed  in 
Section  3-2,  transmission  of  error  code  and  synthesis  of  speech 
afterwards . 

As  before,  the  two  interrupts,  namely  the  speech  and  lineside 
interrupts  are  used.  The  former  obtains  its  signal  from  a 12-bit  A/D 
at  a constant  rate  and  then  output  a 12  bit  processed  sample  for  a 
12-bit  D/A.  The  latter  governs  the  transmission  of  data.  At  every 
interrupt,  it  outputs  a transmitter  clock  pulse  and  decodes  one  bit  of 
received  information.  At  every  8 interrupts,  it  inputs  a new  eight 
bit  receiver  word.  The  program  transfers  control  to  a wait  loop  after 
it  finishes  performing  either  transmitter  or  receiver  functions.  Since 
there  is  no  framing  necessary  in  the  CVS D algorithm,  synchronization 


is  not  needed. 


Figure  3-21.  Block  Diagram  of  CVS  D-D 
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3.^.1  Walt  Loop 

Figure  3-22  show:'.  .1  Plow  chart,  of  the  wait.  loop.  Its 
function,  as  mentioned  earlier,  la  to  use  up  all  the  time  left  over 
after  processing  operations.  It  also  sequences  the  transmitter  and 
receiver  operations  by  examining  the  transmitter  and  receiver  ready 
flops  and  examines  the  front  panel  to  see  if  initialization  is  desired. 

3.4.2  The  Line  Side  and  Speech  Side  Interrupts 

Figures  3-23  and  3-24  describe  the  software  operations 
when  interrupts  occur.  The  function  of  these  interrupts  are 
discussed  in  Sections  3-14  and  3-15. 

3*4.3  Speech  Analysis  and  Synthesis 

Figures  3-25  and  3-26  describe  the  software  operation  of 
the  transmitter  and  receiver  of  CVSD-B.  It  allows  a selection  of 
transmission  rates  (9.6  kBPS  or  16.0  kBPS),  time  constants  in  gain 
filter  (6.3  msec,  10  msec,  20  msec),  time  constants  in  prediction 
filters  (1  msec,  2 msec,  10  msec,  20  msec)  and  compression  ratios 
(10,  20,  60,  100,  150,  200).  Figure  3-27  and  3-28  show  the  software  operation 
of  Forney  and  Qreshi's  CVSD-C.  Figures  3-29  and  3-30  describe  the  operation  of 
CVSD-D  transmitter  and  receiver  where  one  bit  of  memory  is  used. 

3-5  Operating  Procedures  of  CVSD  Program 

3-5.1  Connecting  the  System 

Connect  the  EDM’s  to  each  other  or  through  a modem  as 
detailed  In  Section  3.2.1. 

3-5.2  Reading  in  the  EDM  Software 

Read  in  the  CVSD  program  as  shown  in  Section  3-2.2. 

Initialization  of  the  EDM 

After  the  program  has  been  executing,  reinitialization  can  be 
achieved  by  simultaneously  holding  switch  register  toggle  switch  15 


Figure  3-23.  Interrupt  Loop  Speech  Side 
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Figure  3-24.  Interrupt  Loop  Line  Side 
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Figure  3-26.  Receiver  for  CVSD-B 


2084-76E 


Figure  3-28.  Receiver  for  CVSD-C 


2073-76E 


In  the  down  position  for  each  EDM  and  then  returning  it  to  the  up 
position. 


3.5.M  Program  Options 

The  CVSD  program  enables  selections  of  two  transmission  rates, 
three  slope  adaptation  logics,  six  compression  ratios,  three  time 
constants  for  gain  filter  and  four  time  constants  for  the  prediction 
filter.  The  EDM  switch  settings  and  their  functions  are  summarized 
as  follows: 

i)  Initialization 

BIT  15  down  * Initialization 

up  * program  execution 

ii)  Transmission  Rates 

BIT  7 down  * 16. 0 kBPS 
up  * 9.6  kBPS 


111)  Slope  adaptation  logics 
a)  BIT  8 and  BIT  9 up 


standard  ’-bit  (CVSD-B) 


BIT  2 


BIT  1 


BIT  0 


Compression  Ratios 


BIT  l\ 


BIT  3 


Gain  Filter  Time  Constant 
6.3  msec 


10.0  msec 


1 


0 


20.0  msec 
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BIT  6 BIT  5 Prediction  Filter  Time  Constant 

0 0 1.0  msec 

0 1 2.0  msec 

1 0 10.0  msec 

1 1 20.0  msec 

b)  BIT  8 down  and  BIT  9 up  = Forney  and  Qreshi's  3-BIT  (CVSD-C) 

c)  BIT  8 up  and  BIT  9 down  = Jayant's  1-bit  (CVSD-D) 
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3.6  I/O  Speech  Filters 

Under  this  contract,  GTE  Sylvania  supplied  ^ identical  I/O 
Speech  Filters,  two  for  each  DCA  PSP.  The  I/O  speech  filters 
are  passive,  low-pass  elliptic  filters  that  exhibit  an  essentially 
flat  passband  characteristic  (0-3200  Hz  at  less  than  0.5  dB 
ripple)  followed  by  a very  sharp  transition  region  occuring  just 
below  their  3 dB  attenuation  frequency.  These  filters,  supplied 
by  ESC  Electronics  Corporation,  Pallisades,  N.J.,  are  very  stable 
and  possess  a very  low  noise  figure  as  a result  of  their  passive 
design.  The  mechanical  and  electrical  specifications  insure 
compatibility  with  the  DCA  PSP’s  and  they  can  be  plugged  into 
card  number  1 of  the  PSP  without  any  modifications.  Figure  3-31 
provides  the  frequency  loss  characteristics  of  the  filter  (ESC 
part  number  ESC143F6l)  whose  input  and  output  impedance  is  5000  ohms. 
These  filters  have  a pin  arrangement  assure  that  will  be  inserted  correctly 
into  the  PC  board.  Care  should  taken,  however,  to  not  damage  the 
filter  pins  when  inserting  or  removing  the  filter  from  the  card. 
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devices  04  «rtoilt«4  inh  Into  * flfn*l  *»ll*»i34 

for  OPCQ  tntt«R*nuiton«  Ttv®  p«K#>»«f  4f»l«n  M*  « f44«l«r( 
tui«orl«fttf4  architecture  aKh  MtMy  pr**rm-  eoMrotla4 
parallelism. 

The  form  of  dfilm  will  provide  on  lapleawntatloa  of  W» 
APCQ  coder  which  ui«i  about  dot  of  tiM»  oroettilM  tim®  mil* 
Ohio.  Our  estimates  Indleoto  that  the  AP’Co  algorithm  will  re- 
quire about  Id  mcc.  to  analyte  speech  and  to  *rwcli 

during  each  2%  msec.  data  frame.  Thus,  thf  processor  hardware 
Idles  about  9 Mtc.  out  of  each  frw».  This  Idle  tine  can  te 
used  for  other  telephone  subscriber  functions*  such  as  encryp- 
tion and  telephone  line  control. 

The  cost  of  parts  Is  also  surprisingly  low.  Today's  cost 
for  parts  In  small  quantities  Is  1700  and  their  power  require- 
ments are  Just  over  02  watts.  Thus,  the  APCQ  coder,  while 
mathematically  complex,  can  be  built  today  using  commercially 
available  LSI  processors  for  a relatively  low  cost. 


The  next  section  will  present  the  architecture  of  the  de- 
sign followed  by  a discussion  of  the  LSI  software  Implementation 
of  the  APCQ  algorithm.  Estimates  of  timing,  memory,  cost  and  power 
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?h*  processing  af  u*  <«!*»«  tun  t>* 

tncr«4**a  »>  Adding  M»Mrr  4««r*c$  of  parallelism.  In 

deer,  ten  *.*,  •«  (UwMNto  he#  *pdt*g  «n  («*«««  «r  *n 

•Stern*!  h4M*V«  n^Xti.iirr  or  M*.  «ifl  S«pw»c  n#rfor«t|ftf# 
throat*  )A«r«9i»4  MN*«r«  eomplvstty  ^ tfcrewgi,  * r^rf*n|cr» 
«wmt  of  the  tu*\f.g  ur\ietMr«,  7h*  otUiMUi  Of  tn«>c 
4*»  tgn*  %r«  **  • *r«oMrliM. 


Th#  basic  LSI  ilfMt  proc#sser  Arehlt#«iwre  for  AfC(i  Jn<- 
ptfMft(»tlon  eoMlnt  of  « Mrro  pmH»«r,  • loop  rroe#**«r« 
o^nory  and  peripheral  Tk»ji  *r*  Inimwotcloi  via 

(mo  main  tuifi  **  »Mtn  in  rir*r«  9-1. 

foeh  of  tn#  too  pmtuon  is  d«»io#4  to  perform  tn* 
(••lit  for  which  it  is  best  iuli«4i  u»#  M«r«  processor  for 
••tup  end  program  control , tn#  loop  proc#*»or  (LF ) for  array 
process tn#  in  l(~blt  or  greater  precision  and  lh#  handling  of 
Interrupts  at  a nigh  rate. 

The  nacro  processor  Is  a mini -computer  on  a chip.  The 
Tesas  Instrument  TWS  9900  Is  chosen  as  this  1 6*b 1 t processor . 

It  Is  designed  to  be  capable  of  running  at  a 3 WHc  micro-cycle 
rate. 

The  16-blt  main  memory  contains  both  the  working  storage 
(both  constants  and  variables)  for  the  macro  processor  and  the 
LP  and  the  program  memory  for  the  macro  processor. 

The  loop  processor  Is  the  fast  arithmetic  and  logic  unit. 
As  shown  In  Figure  k-2,  It  consists  of  an  A Feglster  to  Inter- 
face with  the  main  memory  and  the  powerful  ALU  chips,  which  are 
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the  Advanced  Micro  tte¥tefs  AM.*4^!  %-tli  slice  mlreo  pro» 

etMorf. 

Output  of  the  AMI  chip  ('tin  »m*  Mlrimu  M.« 

output  Register  (0  register)  to  the  a***  i.4«  (f»  M)  or  <Arwi|« 
memory  444c*9#  regutor  (»*  register)  Uw  »4«mj  M M), 
Additional  a Sue  to  0 Bus  l«i  Is  pmHH  to  u*  design 

for  Added  memory  access  fl*IIAUlt|. 

The  micro  program  nurtK^ry  i<4  control  MN««r« 
the  Activities  of  the  IT.  The  control  hardware  lncM#» 
micro  program  control  1*4  itorw*  of  ALU  status,  micro  Jump 
control,  subroutine  nesting  end  priority  Interrupt  handling. 

Both  Interrupting  end  non- Interrupting  I/O  ports  sre 
used.  Each  port  CAn  te  Indlvlduslly  Addressed  end  commended 
vIa  the  A Bus,  either  piecing  dote  onto  the  0 Bus  or  Accept  Inf 
dot a from  It.  Interrupting  porte  Issue  service  requests  to  the 
IP,  which  controls  the  buffering  of  Input  or  output  dels  on  s 
priority  Interrupt  besls.  Bon- Interrupting  ports  con  fee  Ac- 
cessed At  eny  time  by  either  the  macro  processor  or  the  LP. 

All  bus-dr lvtng  points  use  three-state  output  devices 
widely  avail Able  In  curr-nt  technology  at  bipolar  speeds.  All 
devices  In  our  design  are  off-the-shelf  LSI,  MSI  and  SSt 
Schott ky  bipolar  chips. 

The  main  ALU,  the  AM  2901.  Is  believed  to  be  the  most 
powerful  micro  processor  chip  available  at  the  present  time. 

It  meets  all  military  specifications  and  can  function  at  a 
clock  rate  of  6.66  MHz,  or  150  ns  cycle  time. 

One  key  feature  of  the  processor  Is  its  16  word  2-port  RAM. 
The  RAM  serves  as  an  orrehip  register  file.  ?-port  makes  It 


easy  to  «ee«u  too  registers  at  *<■*<*  Mae,  la  aggHloo, 
th*r*  t»  ih*  w •'•yliUf,  «i>teti  !i  Klaarlljr  fur 

mull  ft |» | I»m*  *»»«l  ■ • ( *»  t ! r>«tl  Iwa,  Vt^  lel  I ! •'  I M »<<4 

Moi'k  1 1 »/  v ;«a  I;#  In.  1'1,-mi  % - 

Tb#  signal  flow  of  tn*  pr«e***ur  «m  tn  K|fwr«  u,-  ran 
be  •ipl«ln«4  as  following,  input  *«* * ar«  fog  fr«s  ga'e 
aoaory  to  tho  D bus  I Put  a M)  and  )a<«M  by  is*  A 
Thoy  can  cither  b*  save-g  in  «a»  of  tn*  locations  In 
PA*  or  directly  f*m  iMlltf,  «Mch  goo*  Uk*  artfhm** 1*  or 
logical  ep*ratloM.  V output  fro*  tn*  lUf  can  t«  fog  Pars  to 
tho  processor  for  furtnor  All)  manipulat  ion  or  can  t>*  IttfM 
to  tho  0 rogtstor  or  tho  M roglator.  IKt*  dsta  rut  therefore. 


ho  brought  back  to  4ata 


ry  olo  tno  D Hus  or  uaog  a*  a 


memory  address  pointer  to  accoss  no*  goto  from  mo *no>ry. 

Tho  branching  hardware,  which  la  not  shown  In  figure  *•?, 
looks  at  the  result  of  tho  trevtous  cycle.  All  branches  or*  ***• 
cuted  ot  tho  end  of  tho  cycle  following  tho  Instruction  con- 
taining the  branch.  (This  is  tho  result  of  a tit*  lined  retch). 

All  the  operations  are  controlled  by  tho  micro  r^Tcron  Inrtivrilrif  arts 


synchronised  to  1?.JJ  Cte  single 


clock. 
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This  section  presents  the  transnitler  operations  of  the 
APCQ  algorithm.  All  sections  of  the  transmitter  *n»  anal yB8d  and 
sane  are  shram  as  examples.  Overall  timing  and  memory  estimates 
cf  the  entire  APCQ  implementation  are  supplied,  finally,  the 
cost  and  power  of  the  LSI  processor  based  on  counted  chips  are 
documented. 

h . 3 . 1 APCQ  Implementation 

LSI  design  Implementation  of  APCQ  Algorithm  Is  based  on 


th«  of  »*««  aig atrUim-  *+** 

*ttt«r  operation*  i»r«  Us'e-j  to  ftfirt  i.4. 


I/O  !Merr,H 

|AMt  l»riMltM<  loo 

AND? 

• • MUfc  «*|n 

Mi»M4  waveform  »n4  *>o-rn»*| Jf«1  loo 
«x>4 

PAfoo  fd#rn«i«+^  «iut  i ratio* 

*04  4*4v4M*t t*«llOA 

ITOJICIW  CMffUiMl  Ij 

l*»»l  tl 

Error  signal  calculation 

Error  ilr**}  ihovUmiIm 


Figure  »•*  APCtJ  Tr^smlllpf  Ot-r  >1  « 4tl 


Each  U$k  In  the  ir»M*Ulpr  Is  lapi®**®****  a*  sepsra te 
>oft*«rs  for  execution  In  the  loop  processor  (If),  flany  of 
the  software  tasks  Involve  repealed  calculation*  on  dais  array*. 
These  calculations  are  perforated  In  software  loots.  The  macr* 
processor  will  handle  the  loop  anJ  Interrupt  control.  Theref ore, 
the  overall  timing  and  memory  can  be  estimate*  as  those  required 
In  the  LP , the  macro  processor  and  some  overhead. 

The  use  of  the  macro  processor  Is  presented  concept lonally 

in  Figure  ft-5  along  with  the  LP  operational  frameworr  and 
memory  map.  Within  this  organisation,  the  macro  processor  acts 
as  the  master  computer,  and  the  LP  as  a pro* ratable  peripheral 

processor  for  doing  the  complex  computations  required  by  the 
processing  algorithm. 


f)p 

* ® • 

Tills  technique  of  p»Boar<)«  ys«-  hinges  on  the  use  of  two 
j It,  »« **■■'$' i , • •*'»  .*  » tty  i m • 

wtth  software  riwitr*«»|  «i  .«.-•*«  vti  !■  * w.»p  flop.;".  lie 

It tally,  both  “swap  flops"  are  sK  t>  1 and  the  macro  processor 
ts  ersnl«d  access.  The  LP  simply  spina  In  a loop  waiting  for 
one  of  the  SP*s  to  beeone  jfi. 

The  macro  processor  then  sets  up  for  a loop  processor 
“CALL*  by  placing  subroutine  parameters  into  a swap  area  and 
clearing  the  associate**  SH.  The  LF  (when  finished  with  pro- 
cessing  In  the  other  area.  If  any)  will  then  access  the  swap 
area  and  perform  the  called  routine. 

The  macro  processor  (after  checking  for  access  to  the 
remaining  swap  area)  then  sets  up  for  the  next  LP  CALI,  and 
clears  that  “swap  flop*. 

Note  that  the  macro  processor  always  does  all  the  work 
tt  can,  until  both  swap  areas  ar*>  given  over  to  the  LP  and, 
once  given  access  to  a swap  area,  the  macro  processor  accepts 
data  from  It  (when  appropriate)  before  setting  up  the  next  LP 
CALL. 

The  LP,  In  Its  turn,  simply  obeys  the  CALL  when  given 
access  to  a "swap  area",  completing  the  subroutine  (such  as 
an  ANDF  or  Autocorrelation)  and  then  returning  access  to  the 
macro  processor. 

Two  levels  of  Interrupt  processing  are  Indicated  at 
the  LP,  one  for  1NPUT/00TPU?  service  and  one  to  support  data 
transfer  to  the  monitor  Interface  for  test  and  debugging 
operations.  These  operations  are  assigned  to  the  LP  due 
to  its  speed  in  the  first  case  and  due  to  its  ready  access 
to  fast  registers  In  the  second  case. 
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This  framework  provides  for  full  resource  use  within 
the  confines  of  any  given  organization  of  a processing 
algorithm  with  both  the  macro  processor  and  loop  processor 
working  simultaneously  whenever  possible.  With  this  proto- 
col as  background,  we  proceed  below  to  concentrate  on  the 
LP  subroutines  for  the  APCQ  algorithm. 

Examples 

The  functional  Activity  Diagram  for  the  Loop  Processor 
is  used  to  analyze  the  APCQ  implementation.  All  transmitter 
operations  of  the  APCQ  coder  were  analyzed  and  examples  of 
the  more  time  consuming  tasks  are  provided  on  the  following 
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too. 

a . Autocorrelat  lot,  frunct  lot 

Figure  *-f  Is  the  funet  Ion* ! activity  diagram 
of  the  Inner  loop  of  t»,e  In.  autneorrolat  Ion. 

(Mi  1 1 1 f • 1 1 <**it  1 • tit  1 lull  1 : 1 1 I •*«'.|  ty  l«.;»dl*-w* 
the  multiplier  l»it*»  th*  . *'  * I-'**  ’••*>!  t h«  MtltU 
pi  icartd  Into  h t.  The  product  apt ’ *rs  In  ««l  u. 

The  32-blt  result  la  accumulated  In  double  pr#ol- 
3 Ion  in  location  ht>  and  Rf . CO  represents  the  carry 
out  of  the  low-order  accumulation.  Kj,j  la  used  as 
a loop  control  by  presetting  It  to  l?8. 

The  outer  computation  lo  r,  which  la  not 
shown  Tor  the  sake  of  clarity,  performs  the  func- 
tions of  storing  RJ<t  R* , Increasing  R,  up  to  ««*1- 
mum,  setting  R^,  fi.0,  h , • tc. 

The  outer  loop  can  be  performed  In  1$  machine 
cycles  by  the  LP.  The  Inner  loop  takes  22  cycle*. 

Therefore,  the  total  processing  cycle  for  this  func- 
tion Is: 

22  x 128  x 5 ♦ 15  x 5 ■ 1*155  cycles 

Autocorrelation  will  al  be  normalised  such  that 
Cq*1.  Thus,  four  normalisations  will  require  8 multi- 
plications and  1 division. 

The  total  time  to  calculate  the  autocorrelation 
functions  plus  normalization  is: 

14155  + 5 x 22  ♦ 13?  * 14397  cycles 

b.  AMDF 

Pigure  4-7  diagrams  the  inner  loop  of  the  AMDF. 
is  used  to  increase  the  memory  address  by  three. 

R^  and  R^  are  initialized  at  the  memory  address  index. 

The  results  are  accumulated  in  register  pair  and 
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for  iloMblf  pm*lalwi  ^altlw.  CO  represents 
the  earry  oat  of  the  (over  - a*,  %*•  loop 

2»vv4  2»eeu*al?»t<cNS  t wo  ats«|ul«*  |lffer^nr«s. 

It  must  be  employed  I.*  time**  f «*■>  « < **r  «•[*  AWI»F  i<uiiiU. 
ts  usH  its  4 loot*  control , It  Is  preset  at  12. 

The  Inner  loop  requires  11  eyries  per  Iteration. 

The  outer  loop  for  the  AWf  must  repet  the  memory 
s hlr*- as  pointer,  store  results  etc.;  It  requires 
nine  machine  cycles.  Therefore,  the  LSI  processor 
can  Implement  the  AMPK  In: 

13  x 1?  x 60  ♦ 9 x 60  • 9900  cycles 
e.  The  Re  jced  Waveform 

The  speech  signal  l*  rilt»*r«xl  by  a one-tap 
filter  using  th«  pitch  gain,  u,  as  the  tap  coefficient. 
Figure  ft-H  la  the  functional  activity  diagram  of  the 
inner  loop  of  the  reduced  waveform.  R.,  1s  the  memory 
address  pointer  to  pick  up  the  pitch  gain,  a.  Rj, 

Rjj,  are  also  used  as  pointers  to  pick  up  data  from 

memory  or  store  data  back  Into  memory.  The  inner 
loop  takes  2<J  cycles  to  perform  the  function.  15 
more  cycles  are  necessary  to  Initialize  the  memory 
pointing  register  and  some  other  overhead  operations. 
Therefore,  total  required  cycles  are: 

15  ♦ ?K  x 1?8  « 3087 
d . The  Krror  Signal 

Figure  iJ-9  Is  the  functional  activity  diagram 
of  the  Inner  loop  of  error  signals  calculation.  Ry , 

R10  and  are  used  to  accumulate  the  results.  The 

registers  RQ , Rj  , R ^ , R^  are  used  as  memory  address 
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V = S -oS  w 
n n n-M 


I Rj  - HO  - 1 — MA 


n * 1,  2*  • 128 


Branch  exit 
IF  Rs  - 0 


Figure  4-8.  Reduced  Form  Inner  Loop 
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pointers.  is  used  as  a loop  control. 

The  inner  loop  takes  ??«*•?)  cycles. 

The  outer  loop,  which  pcrfems  the  loop  control, 
register  Initialisation  and  stores  the  results 
back  in  memory  requires  an  additional  1?  cycles. 
Therefore,  the  total  required  cycles  are: 

(22  x k ♦ x 1?8  ♦ 1?  i 1?8  * 1ST** 
e.  PARCOR 

Calculation  coefficients  of  the  PARCOR  coeff- 
icients is  a recursive  process.  For  sake  of 
clarity.  Figure  k-10  shows  only  the  cycled  multi- 
plication part,  which  contains  most  of  the  calcula- 
tions. 

Registers  RQ,  R-,  are  used  as  memory  address 
pointers.  R^  Is  used  as  a loop  counter.  Double 
precision  accumulation  ends  up  In  and  R^.  ha  oh 
term  In  the  calculation  takes  two  multiplications 
and  each  of  them  cycles  1,  2 and  3 times.  In  addi- 
tion, two  separate  multiplications  and  fo-ir  divisions 
are  necessary  to  complete  the  task.  30  overhead 
cycles  are  added  to  guarantee  enough  cycles  for 
memory  address  Indexing,  loop  control  and  storing 
results  back  Into  memory,  etc. 

Therefore,  the  required  cycles  are: 

30  x 3 + 20  x (1  + 2 + 3)  ♦ 20  x (1  ♦ 2 x 3)  ♦ 

132  x U + 20  « 878 


Arithmetic  U»ck  Unit 


Figure  4-10.  Partial  PARCOR  Inner  Loop 
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Timing  and  Memory  Estimates 

The  overall  program  liming  Is  dt-i*lv«-d  free  Figure  *-l. 
Total  required  time  f«r  llK*  transmit t«*r  side  la  calculated 
with  some  of  the  calculations  taken  from  the  examples  shown 
In  the  previous  section.  The  result  l*  shown  In  Table  k-l 
below: 


FUNCTION 

TINE 

(«s) 

Input  Normal lzat Ion 

0.71 

I/O  Interrupt 

0.8k 

AMDP 

l.*9 

a,  pitch  gain 

0.8k 

Reduced  waveform  normalization 

0.89 

Autocorrelation  normalization 

2.16 

PARCOR  quantization  dequant- 

ization 

0.13 

Predictor  coefficients 

0.03 

Error  signal 

2.3* 

Error  signal  quantization 

1.17 

Q.  quantization  level 

0.0k 

Total  time 

10.62 

ms 

TABLE  k-l 

The  Estimates  of  Transmitter  Timing  of  APCQ 


In  the  APCQ  Algorithm,  the  receiver  processing  Is  one- 
third  of  the  complexity  of  the  transmitter.  Therefore,  the 
estimated  timing  Is  one-third  of  that  required  for  the  trans- 
mitter. Also,  allowing  15  percent  overhead  time  for  subrou- 
tine control,  data  handling  and  synchronization,  etc.,  then 
the  overall  processing  of  APCQ  algorithm  Implementation  Is: 
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10. 6,?  * Cl  ♦ 11*)  * (1  ♦ 151)  * 16.?  msee 

Memory  used  In  the  US  1 processor  can  t>e  divide  Into 
micro  and  macro  memories.  Micro  memory  h*vl «g  a word  length 
of  bits  Is  used  for  micro  program  storage. 

Macro  memory  having  a word  length  of  )6  bits,  contains 
data  and  the  macro  program. 

The  amount  or  Micro  program  memory  ( ptiM)  is  conserva- 
tively estimated  as  shown  In  Table  a.?: 


Transmitter 

Receiver 

Input  Normalization 

no 

no 

AMDP 

BO 

a,  pitch  gain 

75 

Division  subroutine 

55 

Square  root  subroutine 

65 

Quantization  subroutine 

60 

40 

Dequantization  subroutines 

BO 

B0 

Reduced  waveform  and  normalization 

65 

65 

Autocorrelation  and  normalization 

60 

Predictor  coefficients 

50 

50 

Error  signal 

60 

60 

Initialization  and  synchronization 

100 

Q,  quantization  level 

60 

PARCOR 

100 

Total 

9*0 

355 

TABLE  4-2 

{ 


Microprogram  Memory  Requirements 


the  iwory  as  simtm  In  Table  *.,%  w«- 

should  »Uom  151  ©vernead  for  memory  for  Input  tin*  end  out  - 
putting  the  digital  data.  Therefore,  total  »*»©ry  is: 

<9*0  ♦ 355)  x <1  ♦ 151)  * 1*90 
or  1*90  x *9  * 7)520  bit#.  This  Is  6 chip#  of  )<*K  ho* 

(IK  x 1©  • with  500  spore  words. 

Kaero  program  I/O  Interrupt#,  control  memory  end  dat a 
buffer  for  the  processor  ere  as  shown  In  Table  *-!. 


•| 


* Memories  are  projected  to  1978  technology. 
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DEVI  Cl 


I/O 

OT_. 

00 

6 

. (Ail 

a.  on 

0? 

) 

. *48 

0.6I 

0? 

14 

1.596 

2.57 

08 

6 

0.456 

1.02 

30 

* 

0.076 

0.68 

3? 

2 

0.070 

0.46 

77 

2 

0.720 

0.46 

13? 

1 

0.3*0 

0.59 

133 

1 

0.250 

0.17 

138 

7 

1.575 

5.56 

148 

1 

0.190 

0.58 

153 

2 

0.450 

1.28 

157 

3 

0.750 

1.92 

163 

3 

1.781 

2.37 

m 

9 

4.61 

6.75 

18? 

1 

O.36 

0.68 

2*1 

8 

5.44 

20.56 

251 

2 

0.550 

1.28 

253 

10 

.34 

5.00 

276 

7 

2.10 

3.78 

330 

1 

0.61 

10.00 

37* 

4 

1.8 

10.84 

2401 

4 

3.7 

120.00 

482 

3 

1.426 

30.00 

9900 

1 

1.000 

60.00 

4K  RAM 

6 

3 

36.00 

16  K ROM 

7 

4.2 

101.50 

12-Bit  A/D 

1 

2.43 

195.00 

8-Bit  D/A 

1 

1.6 

49.00 

TOTAL 

119 

41.748 

watts 

1669.36 

Materials , 

Quantities , 

Power  and  Cost 

TABLE  4-4 
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C oat  and  Power  Estimates 

The  cost  and  power  »*at  Iasi*-  of  Hu  |.**|  processor  was  per- 
formed by  adding  u)  th«  •*.>„;»  \ « r n<  . i..  of  >*ueh  elt'i*.  Losses 

caused  by  the  PC  boards,  mechanic  assembly,  power  supply, 
etc.,  were  not  Included  In  the  estimate,  but  those  are  esti- 
mated to  be  small.  The  parte  list,  total  power  and  cost  are 
shown  In  Table  4-4. 

Other  LSI  Processor  Structures 

The  parallelism  of  the  proposed  LSI  processor  can  be 
further  Increased  by  changing  the  loop  processor  as  shown 
In  Ptgure  4-11  to  Include  an  Index  processor  and  hardware 
multiplier.  Now  the  loop  processor  presents  three  func- 
tionally Independent  processors. 

The  Index  processor  Is  used  for  the  computation  of  memory 
addresses.  It  consists  of  a 16-element  file  of  data  plus  an 
8-element  array  of  pointers  with  a United  ALU. 

The  hardware  multiplier  performs  16-blt  by  16-blt  two’s 
complement  multiplication  providing  32-blt  results  In  QH  and 
QL  registers  (16  bits  each).  Once  loaded  and  started.  It  will  complete  Its 
product  in  8 cycles  (or  4 or  2 deoendlng  on  design)  while  other  operations 
can  continue  In  parallel  In  both  the  Index  processor  and  main  AUJ. 

The  main  arithmetic/logic  unit  Is  the  same  as  we  men- 
tioned before.  However,  its  memory  addressing  function  can 
be  performed  by  the  index  processor  and  Its  multiplication 
function  can  be  performed  by  the  multiplier  processor. 

These  three  sections  of  the  loop  processor  run  synchron- 
ously under  microprogram  control. 

In  this  section,  we  consider  adding  either  tne  Index 
processor  or  the  multiplier  processor  or  both,  to  the  main 


ALU.  The  timing,  memory,  cost  and  power  will  be  re-established 
to  compare  with  the  prevlnun  single  main  AUI  processor. 

Timing  and  Memory  Estimates 


With  an  Index  processor,  the  memory  address  pointer  can 
access  memory  locations  at  the  same  time  the  main  ALU  performs 
other  operations.  It  can,  therefore,  save  one  micro  program 
instruction  at  the  beginning  of  every  memory  Index.  With  a 
hardware  multiplier,  a multiplication  requires  10  cycles  with 
an  Index  processor  and  11  cycles  without  an  Index  processor.* 
Without  a hardware  multiplier,  multiplication  requires  20 
cycles  with  an  Index  processor  and  22  cycles  without  an  Index 
processor.  Hence,  the  use  of  a hardware  nultlpller  saves  at  least  half 
the  micro  Instructions.  The  added  parallelism  has  Increased  the  complexity 
of  bus  structures  somewhat  and  the  word  length  of  micro  Instruction  by 
three  bits  for  the  multiplier  and  13  bits  for  the  Index  processor. 

The  timing  of  the  transmitter  side  of  the  APCQ  Implemen- 
tation for  the  different  structures  Is  analyzed  and  listed 
in  Table  4-6. 

The  receiver  processing  time  is  considered  to  be  one- 
third  the  required  time  for  the  transmitter.  Overhead  time 
Is  15  percent  of  the  total  time  required  for  the  transmitter 
and  the  receiver. 

From  Table  4-5,  we  can  conclude  that  adding  a hardware 
multiplier  improves  the  processing  speed  by  about  30  to  50*. 
Adding  an  index  processor  improves  the  processor  speed  by  16 
to  25*. 

* It  is  almost  always  possible  to  perform  additional  useful 
operations  in  either  the  ALU  or  the  index  processor  while 
the  multiplier  is  performing  its  function  such  that  the 
multiply  time  is  effectively  reduced. 


With  InJex  Processors  [ Without  Index  Processors 

With  Multiplier 1 Without  ttiltlpller  [ With  Hiltlpllerl  Without  Hiltl 


TABLE  4-5 

Processor  Time  Requirements  vs.  Hardware  Complexity 
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With  the  different  structure s micro  program  memory  Is 
different  as  la  the  micro  Instruction  word  length  ami  memory 
size. 

Table  — 6 shows  these  differences.  Overhead  cycles 
have  taken  15*  from  the  total  transmitter  and  receiver 
memory. 

Prom  the  table,  It  Is  clear  that  with  the  Indicated 
Index  processor  added,  the  overall  memory  required  is  In- 
creased, because  of  the  Increase  of  micro  Instruction  word 
length.  Other  Index  processors  with  greater  computational 
capability  can  provide  for  greater  speed  Increases  with  ap- 
proximately the  same  increase  in  micro  control  memory. 

The  total  memory,  which  Includes  both  macro  and  micro 
memory  will  be  a3  shown  in  Table  ^-7. 


i 


MEMORY 

With  Index  Processor  Without  Index  Processor 

1th  Multiplier!  Without  Multiplier  With  Multiplier  I Without  Multiplier 
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Estimates  of  Cost  and  Power  of  Index  Processor  r>nd  Multiplier 

The  cost  and  power  required  Tor  the  Index  processor  Is 
$85  and  5.11  watts  respectively.  Kor  the  multiplier  they 
are  $50  and  2.5  watts  respectively. 

Therefore,  the  cost  and  power  for  the  modified  processor 
appear  In  Table  4-8. 

Discussions 

Table  4-5  and  Table  4-8  have  shown  us  many  significant 
results.  The  cost  and  power  do  not  differ  greatly  for  the 
different  structures,  but  the  speed  does.  The  best  choice 
for  APCQ  seems  to  be  the  processor  with  multiplier,  but  with- 
out index  processor. 
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SECT TON  5 

RECOMMENDATIONS  AND  CONCLUSIONS 

5.1  Conclusions 

This  contract  has  resulted  In  the  development  of  a high  quality 
16  Kb/s  speech  digitizer  based  on  adaptive  predictive  coding 
principles  which  Is  robust  to  channel  errors  and  relatively  simple 
to  Implement.  Informal  listening  tests  Indicate  that  the  speech 
quality  of  the  3000  Hz  audio  bandwidth  adaptive  predictive  coder  Is 
noticeably  superior  to  Continuous  Variable  Slope  Del tamodulatlon 
(CVSD)  which  we  also  developed  as  a benchmark.  Experimentation 
with  the  CVSD  algorithms  Indicated  they  could  be  Improved  slightly 
by  adjusting  Internal  parameters,  especially  the  compression  ratio 
defined  as  the  maximum  quantizer  level  over  the  minimum  level. 

Although  this  change  lessened  the  slope  overload  distortion, 
making  the  speech  sound  less  muffled.  Its  voice  quality  was  always 
lower  than  that  of  the  Adaptive  Predictive  coder  at  16  Kbps. 

Using  commercially  available  LSI  electronics  the  16  Kbps 
adaptive  predictive  coder  could  be  implemented  using  about  $700  in 
parts  and  would  consume  about  ^2  watts  of  power.  Our  engineering 
estimates  of  the  hardware  Incorporated  the  Advanced  Microprocessor 
Devices  AMD2901  microprocessor  chip  to  develop  a programmable 
multilevel  adaptive  predictive  coder. 

This  approach,  while  more  complex  electronically  than 
that  used  for  CVSD,  would  result  In  a complete  telephone  terminal 
with  built  In  telephone  line  synchronization,  modem  interface  and 
key  interface  at  about  the  same  cost  as  a terminal  employing  a CVSD 
speech  digitizer.  This  occurs  because  the  programmable  approach 


can  tlmeshare  the  microprocessor  hardware  among  the  speech 
digitizer  and  these  other  tasks  while  the  terminal  with  the  CVSD 
digitizer  must  have  special  hardware  for  each  task.  Because  they 
must  have  this  specialized  hardware,  present  CVSD  terminals  are 
surprisingly  costly,  even  though  the  CVSD  digitizer  Is  extremely 
simple.  Thus,  for  about  the  3ame  terminal  cost,  the  multilevel 
adaptive  predictive  coders  developed  under  this  contract  would 
noticeably  improve  the  voice  quality  over  the  CVSD.  Consequently, 
these  coders  should  be  studied  and  refined  still  further. 

5.2  Recommendations 

The  final  multilevel  adaptive  predictive  coder  system  determined 
under  this  contract  had  a 3100  Hz  bandwidth  and  a 5 level  flxed/frame 
quantizer.  Furthermore,  It  employed  both  a pitch  loop  predictor 
and  a coefficient  predictor.  The  fixed/frame  quantizer  was  used 
because,  at  the  design  review.  It  was  the  only  quantizer  known  that 
was  robust  to  channel  errors.  The  5 level  quantizer  was  used 
because  it  was  the  most  straightforward  approach  to  obtaining  a 
3000  Hz  audio  bandwidth  and  the  pitch  loop  was  necessary  because 
the  speech  quality  degraded  when  it  was  removed. 

Just  as  the  contract  ended,  we  started  to  investigate  modified 
Jayant  quantizers  and  other  quantizers  which  are  also  robust  to 
channel  errors.  Preliminary  investigations  indicate  that  these 
quantizers  provide  improved  quality  over  the  fixed/frame  quantizer, 
and  their  use  can  eliminate  the  need  for  the  pitch  loop  without 
loss  of  quality.  We  recommend  that  these  robust  quantizers  be  studied 
further,  and  that  thev  be  incorporated  into  PSP  software  for  operation 
in  the  DCA  equipment.  This  will  provide  high  quality  adaptive 
predictive  coders  using  these  quantizers.  In  addition,  we  recommend 


that  variable  length  coding  coupled  with  storage  buffers  be  used 
to  obtain  a 3000  Hz  audio  bandwidth  with  either  a 7 or  8 level 
quantizer.  This  will  Improve  the  slgnal>to->nclse  ratio  of  the 
system  while  keeping  the  trar. amission  rate  at  16  Kb/s.  These 
studies  should  simplify  the  Implementation  hardware  from  the  system 
being  delivered  under  this  contract,  while  Improving  voire  quality 
still  further. 
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APPENDIX  A 

MORE  DETAILED  PLOW  CHAHTS  OP  APCQ  CODER  SIMULATION 


A1 


A.l  INTRODUCTION 

Figures  A-l  and  A-2  contain  more  detailed  Plow  charts  of  the 
fixed-point  simulation  with  8 level  Jayant  Quantizer.  These  characters 
Indicate  the  scaling  needed  to  process  the  speech  with  16-blt  arithmetic 
This  scaling  Is  necessary  since  the  summing  of  the  data  can  cause  the 
answer  to  overflow  and  exceed  the  maximum  number  of  the  machine  (2 
Additionally,  scaling  Is  needed  to  avoid  underflow  In  multiplication 
since  only  16  bits  of  a 32-bit  product  of  two  16-blt  numbers  can  be 
retained.  It  Is  common  practice  to  retain  the  upper  or  most  significant 
16  bits,  which  Is  equivalent  to  dividing  the  32-blt  product  by  2J^. 

This  operation  Is  performed  automatically  by  the  hardware  and  Is  called 
a fractional  multiply.  Here  the  numbers  are  thought  of  as  fractions 
having  a magnitude  less  than  1,  and  the  location  of  the  binary  point  Is 
not  changed  after  multiplication.  The  flow  charts  Indicate  this  frac- 
tional multiply  by  showing  the  division  by  21^,  which,  as  previously 
stated.  Is  performed  automatically  by  the  hardware. 

A. 2 MAIN  PROGRAM  PLOW  CHARTS  AND  LISTINGS 

Figure  A-l  presents  a flow  chart  of  the  analyzer  section  of  the 
APCQ  Coder.  In  this  portion  of  the  program,  the  data  is  read  off  the 
tape,  pre-emphaslzed,  normalized,  and  processed  to  compute  the  pitch  M, 
pitch  gain  ALPHA.  This  portion  of  the  program  also  Includes  computa- 
tions and  quantizations  of  the  residual  error  signal  and  the  encoding 
of  the  quantized  error  signal.  Additionally,  the  FARCOR  parameters  are 
calculated  by  this  code. 

Figure  A-2  similarly  describes  the  receiver  portion  of  the 
code.  Here  the  transmission  parameters  are  decoded,  the  predictor 
coefficients  recomputed  from  the  PARCOR  coefficient  and  the  error 
signal  used  to  excite  the  reconstruction  filter.  Prior  to  output  on 
tape,  the  speech  is  de-normallzed  and  de-emphaslzed . 


Figure  A-l.  Flow  Chart  of  Fixed-Point  APCQ  Coder -Analyzer  (Sheet  2 of  3) 
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Figure  A-2.  Flow  Chart  of  Fixed-Point  APCQ  Synthesizer 
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Table  A-l  contains  a lint  of  the  variable  names  and  their  mean- 
ings as  used  In  the  main  program  of  the  APCQ  Coder.  Included  in  this 
Table  Is  a list  of  the  subroutines  and  arithmetic  functions  used  a3  well 
as  the  data  arrays  needed  to  filter  the  speech  data. 

The  FORTRAN  listing  of  the  main  program  is  contained  In  Figure 
A-3  and  Includes  comment  statements  Indicating  the  main  functions  of  the 
program. 


TABLE  A-l.  VARIABLE  NAMES 
SUBROUTINES 

TAPEl(I) 


SW1(I,  ITOG) 

DISREG  (I) 

N0RML1  (V,W,L,N,NS) 

RANDIN  ( I , J , Y ) 

IALPQ  ( IALPHA ) 

IALPD  (IALPHA) 
EDQNT3(MIE,IE(J) , LEVLT1 ) 

RANERR  (A1,A2,A3,A4,A5,A6) 

EDDNT3(IEREC(J) , IER,  LEVLR1 ) 

IPARQ(IPARCT(J) ) 
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USED  IN  MAIN  PROGRAM  OF  APCQ  CODER 
MEANING 

Magnetic  tape  handling  where  I Indicates 
function : 

1=1  read 
1*2  write 
1=3  rewind  tapes 
1=4  search  for  last  file 
1*5  write  end  of  file 
1=6  rewind  only  Input-search  for 
desired  Input  file 

Set  IT0G=1  If  console  switch  I Is  up, 
ITOG*0  if  down 

Display  I on  console 

Normalize  N numbers  in  array  V to  N+l 
bits  & place  result  in  W.  NS=#  shifts 
to  do  this 

Random  number  generator 
Quantizes  pitch  gain  to  3 bits 
Dequant ize  IALPHA 

Code  error  vector  MIE  into  transmission 
vector  IE(J) 

Inserts  random  errors  into  data  stream 
A1  according  to  A2-A6 

Dequantize  error  IEREC(J)  into  IER, 
adjust  Level  LEVLR1 

Quantize  PARCOR  (Reflection  coefficient 
parameter  IPAR  CT(J) 

Dequantize  PARCOR  transmission 
parameter 


j 


IPAROD( IPARCT ( J ) ) 
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TABLE  A-l . VARIABLE  NAMES  USED  IN  MAIN  PROGRAM  OP  APCQ  CODER  (Cont.) 
SUBROUTINES  MEANING 


FUNCTIONS 

IADD  (11,12, IOVR ) 
IFMUL  (11,12) 

IFDIV  (11,12) 

IIMl'L  (II, 12, IOVR) 

VARIABLES 
LEVLTI,  LEVLR1 

LP 

IR1N1  IY1 

NTOTI 

NTUPS 

NTOTO 

NSKIP 

1ST 

NSHFT,  NRSHFT 
LSHFT,  LSHRC 
NIN 
F 

ISPIN 

IV 

IRIN2 , IY2 

IA 

IR 

M 

IALPHA , IRALPH 
MALPHA 


Add  II  to  12,  Increment  IOVR  If  overflow 

Fractionally  multiply  II,  12;  round 
result 

Fractionally  divide  II  by  12 

Integer  multiply  II  by  12;  Increment 
IOVR  If  overflow 


Error  quantizer  level  at  transmitter 
and  receiver 

Maximum  value  of  quantizer 

Coefficient  predictor  filter  memory  and 
transmitter  receiver 

Tape  handler  data 

Tape  handler  data 

Tape  handler  data 

Tape  handler  data 

Tape  handler  data 

Present  normalizer  shift 

Past  frame  normalizer  shift 

Speech  Input  array 

Input  array  after  pre-emphasis 

Input  array  after  normalization 

Reduced  waveform 

Pitch  filter  history  (reconstructed 
waveform) 

Predictor  coefficients 
Autocorrelation  values 
Pitch  delay 
Pitch  Gain  Alpha 
Quantized  Alpha 
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TABLE  A-l.  VARIABLE  NAMES  USED  IN  MAIN  PROGRAM  OF  THE  APCQ  CODER  (Cont.) 
VARIABLES  MEANING 


MIE,  MRE 
IPARCT,  IPARCR 
ITRAN,  IEREC 
YY 
Y 

NOUT 


Quantized  error  signal  array 
PARCOR  parameters 
Coded  vector  of  error  signal 
Output  array  before  de-emphasis 
Output  array  after  de-empnasis 
Output  array  to  tape 


A-2 . QUANTIZATION  OF  TRANSMISSION  PARAMETERS 

Table  A-2  shows  the  quantization  of  the  filter  parameters,  the 
normalization,  and  the  error  signal. 


TABLE  A-2.  NUMBER  OF  TRANSMISSION  BITS 
FOR  EACH  PARAMETER 


Parameter 

No.  of  Bits 

pitch  gain 
pitch 

normalization 

PARCOR 

error 

3/frame 
6/frame 
3/frame 
12/frame 
2 . 1/sample 

The  quantization  and  dequantization  of  the  pitch  gain  ALPHA  is 
performed  by  the  subroutines  IALPQ  and  IALPD  respectively.  Quantiza- 
tion and  dequantization  of  the  PARCOR  coefficients  is  performed  by 
subroutines  IPARQ  and  IPARDD  respectively,  with  each  PARCOR  quantized 
to  ^ bits.  Quantization  of  the  pitch  is  performed  simply  by  subtracting 
15  from  the  calculated  value  and  quantization  of  the  normalization  by 
subtracting  1 from  the  unquantized  value.  The  quantized  values  then 
run  from  0-60  for  pitch  (6  bits)  and  0-7(3  bits)  for  normalization. 
These  numbers  represent  coded  values  for  transmission. 
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A. 3 QUANTIZATION  OK  THE  ERROR  SIGNAL 

Quantization  of  the  error  signal  In  accomplished  by  means  of 
an  eight-level  adaptive  quantizer  as  described  previously  in  Section 
2.8.2.  This  quantizer  technique  was  developed  by  Jayant  as  described 
in  the  IEEE  Proceedings  of  May  197**. 

Figure  A-4  shows  a flow  chart  of  the  subroutine  EDQNT-3  that 
is  used  to  implement  the  adaptive  quantization  process.  This  routine 
converts  the  error  signal  into  one  of  eight  numerical  values,  and 
adjusts  the  quantizer  step  size  in  accordance  with  the  "Jayant 
algorithm" . 

Table  A-3  shows  the  variable  names  and  their  meanings  as  used 

in  subroutine  EDQNT-3.  Figure  A-5  is  a FORTRAN  listing  of  this 

• ) 

subroutine . 


TABLE  A-3.  VARIABLE  NAMES  USED  IN  SUBROUTINE  EDQNT3 


VARIABLE 


MEANING 


LP 

MIE 

LEVEL 


Maximum  value  of  quantizer 
level 

Input  unquantized  error  and 
output  quantized  value 

Present  quantizer  level 


1 
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C 

c 

c 


10 

c 


c 

3 


C 

2 


C 

1 


SUBROUT I NE  EDQNT3 (IN,  I ONT .LEUEL ) 
ERROR  QUANTIZER  AT  THE  TRANSMITTER 
EIGHT  LEUEL  ADAPT  I UE  QUANTIZER 
u i tw  rnFFFiriFNTS  AS  PER  JAYANT 


DIMENSION  C0EFC4) 

DATA  LP/1024/ 

DATA  COEF/0.9,0.9, 1.25. 1.75/ 


1SIGN=1 

IFdN.GE.01  GO  TO  10 

IN--1N 

ISIGN=-1 

IFCIN.LE. LEUEL)  GO  TO  1 
IF(IN.LE.2*LEUEL>  GO  TO  2 
IFC IN.LE.3*LEUEL>  GO  TO  3 
TOP  MOST  QUANTIZER  LEUEL 


1QNT=0 

IFdSIGN.EQ.ll  I QNT = 7 
1N=3.5*LEUEL*ISIGN 
LEUEL-LEUEL*COEF(  -4 ) 

I F ( LEUEL . GT . LP ) LEUEL *LP 
RETURN 

SECOND  FROM  TOP  LEUEL 
j QNT  * 1 

lFdSlGN.EQ.  1)  I QNT =6 
IN»2.5*LEUEL*ISIGN 
LEUEL  « LEUEL*COEF  C 3 1 
IF(LEUEL.GT.LP)  LEUEL=LP 


SECOND  LEUEL  FROM  THE  BOTTOM 
IQNT=2 

IFdSIGN.EQ.l)  IQNT-5 
IN'1.5*LEUEL *1 SIGN 
LEUEL= LEUEL *C0EF(2) 

IF (LEUEL. EQ.0)  LEUEL- 1 


RETURN 

BOTTOM  LEUEL 
I QNT  *3 

IFdSIGN.EQ.ll  IQNT-4 
IN«LEUEL/2XISIGN 
LEUEL=LEUEL*COEF( 1 1 
I F ( LEUEL . EQ . 0 1 LEUEL =1 


RETURN 

END 


Figure  A-5:  FORTRAN  Listing  of  Error  Signal 

Quantizer 


m *'"7 

nu  « 

\ liUUuU  i • 


Ji 
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A. 4 Dequantization  of  the  Error  Signal 

The  flow  chart  of  Figure  A-6  shows  the  dequantization  of  the 
error  signal  via  subroutine  EDDNT3.  This  algorithm  performs  the 
inverse  operations  to  those  performed  by  the  quantizer  at  the  analyzer 
and  yields  as  output  the  sequence  of  quantizer  error  levels  generated 
in  the  analyzer. 

Table  A-4  contains  the  variable  names  and  their  meanings  used 
in  this  subroutine.  A FORTRAN  listing  of  this  subroutine  is  contained 
in  Figure  A-7. 

TABLE  A-4 . V/RIABLE  NAMES  USED  IN  SUBROUTINE  EDDNT3 
VARIABLE  MEANING 

LP  Maximum  value  of  quantizer  level 

MIE  Input  quantized  value 

ISIG  Output  dequantized  value 

A . 5 Other  Subroutines 

A. 5-1  Introduction 

This  section  presents  various  subroutines  that  have  not  been 
discussed  previously.  Table  A-5  presents  a list  of  the  names  of  these 
subroutines  and  their  functions. 


. - 
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DEQUANTIZE 
ERROR  SIGNAL 
EDDNT3 


; 

ISIGN  - -1 

UN  =■ 

IN+  1 

S' 

NO 

2061-76E 


Figure  A-0.  Flowchart  of  Subroutine  ED0NT3  - Dequantizing  Error 
Signal  (8-Level  Jayant  Quantizer) 
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SUBROUT 1 1 IE  £DDNT3( I N . 1 01 JT , LEUEI . ) 

SUBROUTINE  FOR  DEQUANT I 2 I NO  3 BIT  ADOPT I UE 
QUANTIZER  SIGNAL  AT  THE  RECEIVER 
FULLOUS  THE  JAYANT  ALGORITHM 


10 

C 

1 

11 


2 

3 

33 

4 


DIMENSION  COEF < 4 1 

DATA  COEF/ 1.75.1.25.0.9,0.9/ 

DATA  LP/ 1024/ 

ISIGN—1 
I IN* IN+1 

I F ( I IN. LT. 1 . OR. I IN.GT.8)  I IN=1 
IF< I IN. LE. 4)  GO  TO  10 
ISIGN=1 
I IN=9~I IN 

GO  TO  (1.2,3, 4) I IN 


I OUT-3.  SKLEUEL* (SIGN 
LEUEL=LEUEL*COEF( I IN) 
IF(LEUEL.GT.LP)  LEUEL*LP 
RETURN 

I OUT *2 . S*LEUEL* I S I GN 
GO  TO  11 

I OUT1 1 , 5YLEUEL* 1 5 1 GN 
LEUEL * LEUELfcCOEF ( I IN) 

I F ( LEUEL . EQ . 0 ) LEUEL= 1 
RETURN 

I OUT-0 . 5*LEUEL *1 5 I GN 

GO  TO  33 

END 


Figure  A-7 : FORTRAN  Listing  of  Error 

Signal  Dequantizer 


c ' 


TABLE  A-5. 
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SUBROUTINES  AND  THEIR  FUNCTIONS 

SUBROUTINE 

NAME 

FUNCTION 

IALPQ 

Quantize  Alpha  To  3 Bits 

IALPD 

Dequantize  Alpha 

IPARQ 

Quantize  PARCOR  coefficient 

IPARDD 

Dequantize  PARCOR  coefficient 

RANERR 

Generate  and  Introduce  Random 

Errors 

N0RML1 

Normalize  Input  Speech 

IADD 

Fixed-Point  Addition  With  Overflow 

Indication 

ISUB 

Fixed-Point  Subtraction  With 

Overflow  Indication 

IFMUL 

Fixed-Point  Fraction  Multiply 

With  Rounding 

ITMUL 

Fixed-Point  Integer  Multiply  With 

Overflow  Indication 

IFDIV 

Fixed-Point  Fractional  Divide 

5.2  Subroutines  IALPQ 

and  IALPD 

Subroutines  IALPQ  and  IALPD  each  contain  one  argument.  In 
IALPQ  the  argument  upon  entry  is  the  unquantized  ALPHA.  The  sub- 
routine converts  this  argument  to  a three-bit  representation  having 
values  0 through  J.  These  values  are  in  a form  suitable  for  trans- 
mission. Subroutine  IALPD  accepts  the  corresponding  received  three- 
bit  value  and  converts  it  into  one  of  eight  possible  values  for 
ALPHA  by  a table  look-up  procedure. 

Figure  A-8  contains  a FORTRAN  listing  of  both  IALPQ  and 


IALPD. 


sSKi™Rbx?tf[j),NG  "LPH0  ™D  DKQunKr,z,NS 


IF'J.LT.0)  J—J 
IFlJ.GT. 23250)60  TO  1 
1 F( J. ST. 31550) GOTO  2 
IF(J.ST. 30650)  GO  TO  3 
lFfJ.GT.23PS0)GO  TO  A 
1F<J.GT.2S8&?>  GO  TO  5 
I F(J. ST. 20*100)  GO  TO  & 

IFfJ.GT.il 850  > GO  TO  7 

J-7 

RETURN 

J-0 

RETURN 

J = 1 

FETURI  1 

J-2 

RETURN 

J-3 

RETURN 

J**4 

RETURN 

J*5 

RETURN 

J-6 

RETUPN  ■ , 

END  ! 

PROGF.NM  FOP  DEQUANT 1 2 1 NG  ALTHfi 
SUBROUTINE  iHLPDf Ji 
DIMENSION  lOLPfS) 

luj+l  ' r<LP  3??S0’ ?17r'1  • 309S1 . 29701 . 27301 , 2297b.  15978.5850^ 

J-IALPfJ) 

RETURN 

END 


Ciipy  AVAIL/iOLE  T9  BDC  ElLSJ'ST 

PERMIT  FliLLl  LEGisLL  fitCuoUiiUu 


Figure  A— 8 : FORTRAN  Listing  of  Subroutines  IALPQ  and  IALPD 
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A. 5- 3 Subroutine  IPARQ  and  IPARDD 

Subroutines  IPARQ  and  IPARDD  quantize  and  dequantize  the  PARCOR 
coefficients.  Each  subroutine  takes  one  agreement  and  returns  the 
answer  in  this  same  variable.  In  IPARQ  the  argument  upon  entry  is 

i. 

the  i unquantized  PARCOR(l).  The  subroutine  converts  this  in  a 
A bit  representation  having  values  0 through  15.  These  values  are 
In  a form  suitable  for  transmission.  Subroutine  IPARDD  accepts  the 
corresponding  4 bit  value  and  converts  it  into  one  of  16  possible 
values  for  the  ith  PARCOR  coefficient.  Figure  A-9  contains  a FORTRAN 
listing  of  both  IPARQ  and  IPARDD. 
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SUBROI  ITI NE 

I FLN6=0 
IFf J.GE.0) 
J=lHBS< J) 
J= J+16 
RETURN 
END 


IPARGKJ) 
GO  TO  1 


SUBROUTINE  IPHPDD(-J) 

DINENSION  I Pin' 16) 

DOTH  IPD'1O00. 3100. 5075 » 7025 . 9325 » 11200,13250, 

15225, 17500. 19475, 21450, 23400, 25700, 27075. 29550, 31500X 

J s J ' 1 

rFLAG=0 

IF( J.LE. 16)  GO  TO  1 

IFLAG=1 

J= J-16 

J'IPD(J) 

IF( IFLftG. EQ. 1 ) J=-J 

RETURN 

END 


Figure  A— 9 : FORTRAN  Listing  of  PARCOR  Quantizer 

(IPARQ)  and  of  PARCOR  Dequantizer 
( IPARDD ) 
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A . 5 - ^ Subroutines  RANERR  and  IXOR 

Subroutine  RANERR  is  entered  with  six  arguments  for  introducing 
channel  errors  into  the  transmission  parameters.  These  arguments 
include  the  transmission  data,  the  probability  of  channel  bit  error, 
states  for  the  random  number  generator,  and  a count  of  the  total  number 
of  errors.  The  assembly  language  program  IXOR  is  used  in  conjunction 
with  RANEPR  to  exclusively  "OR"  the  transmission  data  with  the  random 
bit  error  data  generator  within  RANERR.  Documentation  of  subroutines 
RANERR  and  IXOR  are  contained  in  Table  A-6  and  Figure  A-10.  Table  A-6 
lists  the  variable  names  contained  in  subroutine  RANERR  and  their 
meanings  as  used  in  this  program.  Figure  A-10  contains  listing  of  the 
FORTRAN  subroutine  RANERR  and  the  assembly  language  subroutine  IXOR. 

TABLE  A-6.  VARIABLE  NAMES  USED  IN  SUBROUTINE  RANERR 


VARIABLE 

MEANING 

MASK 

Mask  for  Introducing  a bit  of 
error  in  either  first  through 
8 bit  of  word 

RANDU 

Subroutine  for  obtaining  random 
number  YRN  between  0 and  1 

PROB 

Probability  of  error  any  bit 

IXOR 

Function  for  exclusive  OR  in  bit 
of  error 

I WORD 

Output  word  with  errors 
introduced 

SUBROUTINE  RLMERR1  I UOPD.  KEYM.  PROB.  IPN.  JRN.KCTJNT) 

C IUORD  - THE  UOPD  CONTAINING  THE  CAJANT I TED  PARAMETER 
C KEYM  - WHICH  BIT  OF  IUORD  THE  SUBROUTINE  WILL  OPERATE  ON 
C PROB  - THE  FRACTION  OF  ERROPS 

C JRN  PND  JRN  - NOS.  USED  Bi'  THE  RANDOM  NO.  GENERATOR 
C HUNT  - A RLWING  TOTAL  OF  THE  ERRORS  IN  THIS  FRAME 


C 

C A RANDOM  ICMBER  BETWEEN  0 AND  1 IS  GENERATED  AND  IF  IT 
C IS  LESS  TWN  PROB.  IWORD  IS  EXCLUSIUE  OR-ED  WITH  THE 
C APPROPRIATE  MASK  IN  ORDER  TO  FLIP  THE  BIT  POINTED  TO  BY  KEYM. 


C RANERR  IS  CALLED  OTCE  FOR  EACH  SIGNIFICANT  BIT  IN  IWORD. 
C I.E.  IF  A PARAMKTER  IS  QLWNTIZEP  TO  4 BITS.  RANFRR 
C IS  CALLED  WITH  KEYM  EQUAL  TJ  1.  TO  2.  TO  3.  AND  TO  4. 


DIMENSION  MASKI0) 

DATA  fVSK  .'1.2.4.8.16.32.64.126/ 
CALL  RANDLUIRN. JPN.YFtl) 

IF  (YRN-PPOB)  2O0.SOO.SOO 
200  JUNK  - I XORi IWORD.  MASK (KEYM;) 

I LORD- JUNK 
KOUNT’KOUNT+1 
900  RETURN 

END 


.TITLE  IXOR 
. GLOBL  IXOR 

I THIS  ASSEMBLY  LANGUAGE  SUBROUTINE  TAKES  EXCLUSIUE  OR  OF 

ITWO  INKJT  ARGUEMENTS  AND  LEAUES  hNSWEP  IN  REGISTER  0 FOR  A FORTRAN 

SFUCTIOH  CALL 

R0-V.O 

Rl-Yl 

RS-P-S 

'iXOR*  TST  (RS)+  »Jlf1P  OUEP  B RAICH 

rcu  0 iRS ) + . RO  i GET  FIRST  APGUEMENT 

MOU  0(RS)  + .R1  tGET  SECOND  ARGUEI1ENT 

XOR  Rl.RO  i EXCLUSIUE  OR  U1TH  ANSWER  IN  RO 

RTS  RS  iLEAUE 

• END 


I 

I 


Figure  A-10:  FORTRAN  Listing  of  Subroutine  RANERR  and 

Assembly  Listing  of  IXOR 


A. 5-5  Subroutine  N0RML1 


w ■ 


i 

| 


I 


Subroutine  N0RML1  contains  five  arguments,  the  unnormalized 
Input  array,  the  normalized  output  array,  the  n .mber  of  data  points 
In  the  array,  the  value  of  the  largest  point  af-  n rr.a . 1 zat  Ion , and 
the  number  of  shifts  computed  by  the  subroutln*  . i.*»;s*mbly  language 

listing  of  this  subroutine  Is  contained  In  Pl,-.r< 

A. 5. 6 Function  IADD , ISUB,  IFMUL,  IIMUL  and  IFDIV 

The  routines  discussed  In  this  section  are  called  "Func* Ions" 
since  they  return  an  argument  to  the  main  program  via  a specified 
register  rather  than  one  of  the  calling  arguments,  functions  IADD  and 
ISUB  arguments  each  contain  three.  The  first  two  are  the  Input  numbers 
to  be  added  (or  subtracted)  while  the  third  argument  Is  incremented  If 
an  overflow  occurs.  The  function  IFMUL  has  two  arguments  and  returns 
the  upper  sixteen  bits  of  the  product  (with  rounding)  of  these  two 
arguments  to  the  calling  program.  This  function  is  sometimes  referred 
to  as  a fractional  multiply  because  it  preserves  the  location  of  the 
binary  point. 

The  function  IIMUL  has  three  arguments  which  are  the  two  input 
numbers  to  be  mrltiplied  and  the  overflow  indicator.  This  function 
returns  the  lower  16  bits  of  the  32-bit  product  and  is  referred  to  as 
an  integer  multiply. 

The  function  IFDIV  accepts  two  inputs  and  returns  the  frac- 
tional part  of  the  division  to  the  calling  program. 

Assembly  listing  of  these  fixed  point  arithmetic  functions  are 


given  in  Figure  A-12. 


. TITLE  Nnwt.X.P«L 

SOUTINE  r-firn  t.tc  < r.  shifting  to  the  left  >>x  hfppy  of  integer 

hUTHEP-:  OM  *|I«D  PP-H  H FijETPfH  HpRPY  Mll>  PLACES  THE  RESULT  IN  p 
HEU  PKRP.  UITh  .jiE  u)FD  PER  POINT  INSTEAD  OF  THE  TUT  UQPDS 
►XiaiEP  IN  FUFTRW) 

CY4.L  IOPIL' IN.LUJT-N.NSPEC.riS' 

IH-FOPTPMI  INIUT  UP PPY 

CUT •«'<  n PUT  NPPP < • 0»IE ' IIJRIV  INTEGER ) 

14*  « "F  UW  IN  Apr", 

N.FE'  *SPE>  IFICmTIiJI  l'P  SHIFTED  OUTPUT 
NS*»  UF  SHIFT  CIjnrvjTEU  III  f«iFM 
.GLUCL  NUFTV  1 

Pi>  • *> 

K 1 I 
K2*  jr 
p>  • i 

P-4*-  -I 
PS*  s 


I 


FOFTV  1 i 

TST  '£•?>♦ 

i SET  OUER  FOE TEN  PPM  O' 

IOU  fM.FVS 

MOU  Fi.Pis. 
n:n  p.?.s 
noo  P3.P3S 
MO"  F4.F4S 

iCAUE  GENEFVX  REG  1 STEP'S  0. 1 .2.3. -»  fcS 

'■rj.i  . P5'4.po 

iGET  PDPE9S  OF  INPUT 

IT?.1  SO.  I NSW 

ibPUE  INPUT  hDPESS 

MC"  ' FS  • PI 

iOlTPUT  gdfess 

fEi"  *•  FS>*.P2 

: a OF  MTP  UOKDS 

MC *»  P2.H 

: Su>  IE  ■ OF  U3FDS  IN  PRRPY 

FftU  •f'FS’+.FO 
'O'  P3.NS.r-EC 

iSPUE  LOCATION  OF  SHIFTED  UPLUE 

CLP  P-J 

iJNCLUSI'.'E  OP 

LOOP  • 

N»  i • PC  • ♦ . P * 
PFL  Fd/i 

"SET  1 1 IPLIT  POINT 

NEC-  P3 

'CPH'T  PIS  HE'S.  »-S 

M US' 

PIS  po.p4 

i P4  CONTAINS  SUM  OF  It-Cl.USI'/E  OR 

TST  «PO»4 
SOB  P£*LOOP 

iTUO  WORD  IMTEGEP 

CLP  S3 
TST  M 

■ SET  up  to  cam  ■ of  shifts 

PEU  hLI 

' F'PFOEt ITS  CCHTIUJCJUS  L COPING  ONE  PLL  0’S 

LOOM  ' 

lit:  F.< 

l/'-l.  F4 

'SHIFT  LEFT  OI<  E 

BUT  i.iTOPl 

i CHEL I FUP  OUEFFl  OH 

i«  OF  LEFT  SHIFTS  PFQ'J  1 PEP 

ALL' 

IP'  H.R2 

sHOFUHLirE  INPUT 

tIJ  ‘ »JS.  .P" 
SUE  NSFEC. PC 
SUB  F\-.R? 
HOU  IlfW'.Ri' 

i SPECIFY  LCt.  OF  DEC. POINT 

LOOP.?  * 

f 1L»>  J ' PUN.P4 

'GET  INPUT 

T ST  > FV'4 
HSU  P?.P4 

sONB  HOPE  IfHEGEF 

MOU  P-». • p t 1 ♦ 
TST  'Rl'4 
SOP  PS  • LOOPS' 

'DEPOSIT  INTO  OUTPUT 

MO'  • P3  ■ 1 F S • * 
MOU  P.’S • PC 

nou  Fis.  pi 

NUU  FES.  PS 
MC-1  F'Cl-F. 

IO'  * I3.pl 

iretupn  shift  cam 

FT'S  FS 

! PEP  JPN 

i £3  JEFFS' 

STOP.  *-F. 

Pv>i  : 

UIOPD  u 

"SENFPHL  REG  I STEP  TENPORPRY  ST0PP6E 

PIS' 

.HOP'D 

P’S' 

.UOFL  •:> 

P3S' 

.NjpD  0 

P4S' 

.uupr  u 

N: 

.wop:  o 

iLENGH  OF  INPUT  PRRP ' 

IMSPNi 

.nop:  o 

'hDPES.S  Or  INPUT  hSPmi- 

NSFEC  ■ 

. i4jp;  u 
. END 

i -.PEl.  of  SHIFTED  a.TPUT 

Figure  A-ll:  Assembly  Language  Listing  of  N0RML1 


-TITLE  I (WITH. MAC 


i 


Figure 


A3? 


9 

• this  function  adds  two  numbers  pnp  increments  the 
iThikd  argument  ie  there  is  an  ouerflou 

“SAME  FOR  THE  SUITFACT  CASE 


l 


RO*W  1 

Pi*;n 

FS--.S 

-C-4.0PL  I ADD*  I SUB.  IFMUL.  IFDIU.  1 1IUL 
(A* I ADD! I . J. F ' 

: IA-I+J 

■K-K+i  IF  OUERFLOU 


I 


IADD> 


ft: 


TST  '.  RS)  + 

MO)  e'RS)*.R0 
ADD  ft'  PSi  + .RO 
BUS  ft 
TST  (RS)  + 

RTS  FS 
INC.  ft' RSi+ 

PTS  FS 


i GET  ODER  BRANCH 


(ADD 

:JLMP  ON  OUEPFLOU 
iGET  OUER  ARGUMENT 


i 


ISiJB*  TST  'RS) 

MOU  ft'RS)-».R0 
SUB  ft<RS)+.RO 
BUS  A 
TST  (PS)* 

PTS  PS 

(THIS  SUBPOUTINE  MULTIPLIES  TWO  NUMBER  IN  FRACTIONAL  FORM  WITH  ROUNDING 
lAS  IN  THE  GTE  SYLUANI A PSP  INSTRUCTION  MLY 
IFMUL:  TST  i RS)  + 

MOJ  ft ( RS ' ♦ » PO  (GET  THE  FIRST  NUMBER 

MUL  ft(RS)  + .R0  IMULTIPUY  BY  THE  SECOND (33  BIT  ANSWER) 

flSHC  t»l.R0  i EL  I M I NATE  ONE  SIGN  BIT 

TST  R1  » ROUND  IF  THE  17TH  BIT  IS  A ONE 

BGE  NPND  (NO  ROUNDING  NECESSARY 

IMC  R0  iROUNTING  OCCURS  IF  R1  IS  ICGATIUE 

riKTiP  * rTS  PS 

(FUNCTICN  IFDIU: DIUIDEND. DIUISOR) 

IFRACTICMAL  DIU1DE 


IFDIU: 


T3T(RS)+ 

MOU  ftiRS.'  + .Rft 
CLR  R1 
ASHC  K-1.P0 
DIU  ft'  R5)+.R0 
RTS  RS 


i GET  DIUIDEND 


l 


i ADJUST  FOR  DIUIDE 
(DIUIDEND/DIUISQR 
i RETURN 


iFUNCTTCN  IIMUL 
l INTEGER  MULTIPLY 
i IMUL(ft.B.C) 

(ANSWER  - A*B  WITH  C INCREMENTED  IF  OUERFLOW 
I 

TST  (RS)* 

MOU  #(RS>*.R0 
CLC 

MUL  ft(RS>*.P0 
MOU  R1.R0 
BCS  A 
TST  (RS)* 


I IMUL • 


(GET  FIRST  NUMBER 
iCLEHP  THE  CARRY  FLAG 
: MULT I PLY  THE  NUMBERS 
(MO1  INTEGER  PART  TO  V0 
i BRANCH  IF  OUEWLOU 


RTS  RS 
.END 


A-12:  Assembly  Listing  of  Fixecl-I'o! nt  ArHtimeU':  KuncUontJ 


